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ABSTRACT 



Constant load creep data were obtained for four plastics: polyethyl- 

ine, polypropylene, nylon, and polycarbonate. Viscoelastic responses of 
these materials proved to be non-linear in both time and stress* An 
empirical creep equation, based on a time dependent retardation rate 
constant, was developed to characterize the creep behavior and good cor- 
relation was obtained for various stress levels. By assuming the de- 
veloped relationship to be an equation of state, creep behavior of simple 
composites of plastics was predicted with good results when viscoelastic 
responses of the components were similar. The importance of stress-time 
history is demonstrated by using the creep equation to estimate mechanical 
behavior in other modes of testing. The concept of a retardation spectrum 
and its application to the materials studied is considered. 
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NOTATION 



'T - shear stress 

V - shear strain 

- viscosity coefficient 

Y ” shear strain rate 
G - shear modulus 

E - tensile modulus 
t - time 

E - engineering strain 

- strain rate 

CT - engineering stress 

• 

O” - stress rate 
<%- initial stress 
(Tt- stress at time 11 'Ll m 
3to- creep compliance 
0 - current time 

- natural logarithm base 

- glass transition temperature 

- free energy of deformation 

Eo- zero time strain (instantaneous elastic strain) 
delayed elastic strain 
Er(^)“ relaxation modulus 

- non-recoverable flow strain 
S - deformation 

£ - deformation rate 
PE - polyethylene 
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pp - 
PC- 
NY- 
T - 
b - 

flfc') - 
L(r). 

£ - 
A,- 
Az- 
A,- 

C - 

Q- 

K,e>,P- 



polypropylene 

polycarbonate 

nylon 

retardation time constant 
rate constant for delayed strain response 
rate constant distribution function 
retardation time constant distribution function 
delayed strain fraction 

creep parameter representing zero-time strain 

creep parameter representing steady-state creep rate 

creep parameter representing magnitude of delayed 
elastic strain component 

time dependent rate constant in empirical creep 
equation 

time exponent in empirical creep equation 
stress coefficients in generalized creep equation 
stress exponents in generalized creep equation 
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1 . 



Introduction. 



Composite materials with high strength-to-weight ratios show great 
promise. [3] Composites may be looked upon as intimate mechanical combina- 
tions of two or more materials, but differ from multiphase materials such 
as metallic alloys in that they are not homogeneous on a macroscopic scale. 
Composite materials are not new. Concrete and filled resins, for example, 
have been successfully employed for years, and currently a bewildering 
number of combinations of many thousands of materials is available. Of 
necessity, however, the work with composites to date has been directed 
towards development of materials to meet specific technological require- 
ments, and very little research has been carried out in the area of 
characterization of the behavior of composites. [1, 2] 

Physical properties often constitute the starting point for the 
materials selection process, and so characterizing mechanical behavior 
of composites is an important first step in a fundamental research study. 

Of particular interest is viscoelastic behavior, frequently manifested by 
the component materials of the composite structure. Deviations from 
true elasticity are negligible for engineering purposes for materials 
such as metals and glasses, but in polymeric systems, which are extensively 
employed in composites, behavior is dominated by viscoelastic phenomena. 

[2] Many conventional composite materials have mechanical properties 
that exhibit appreciable time dependence through creep and relaxation pro- 
cesses. [4] The stress analysis becomes quite complicated, but it must 
be performed if a viscoelastic composite is to be utilized in an effici- 
ently designed structure giving satisfactory performance in service 
environments. [5] 

In this investigation the main purpose was to characterize the creep 
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behavior of individual viscoelastic materials by analytic equations, and 
then obtain suitable correlations so that the behavior of composites of 
these materials could be predicted. . In addition, these analytic expressions 
were used to estimate performance of materials in modes of testing other 
than creep, such as stress relaxation and constant strain rate uniaxial 
tension tests. 

For the experiment, four different thermoplastics (polyethylene, poly- 
propylene, nylon, and polycarbonate) were chosen, and composites were 
formed as layered sandwich structures. A simple composite was selected 
for this basic study to eliminate as many variables as possible, since 
stress analysis of polymeric materials is inherently complex, even if 
linear viscoelasticity is assumed. In all experiments the specimens were 
subjected to uniaxial tension, and it has been shown [6] for this type of 
loading that changing the type and amount of component interfacial con- 
tact has negligible effect on composite behavior. This investigation was 
carried out in a laboratory with controlled temperature of 23 + 1° C. 
and relative humidity 50 + 2% in accordance with ASTM Standards. [7] 
Laboratory set-up and equipment is reported in Appendix G. 
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2 . 



Linear Viscoelasticity. 



For the ideal elastic solid, stress is proportional to strain when 
a load is applied in accordance with Hooke's Law. It is also independent 
of the rate of strain. For example, at stresses below the yield strength, 
the behavior of most metals is given by Hooke's Law: 



(2.D 7~= gy 



where = applied shear stress 
'T = shear strain 
G = modulus of rigidity 

For the ideal viscous liquid, stress is proportional to the rate of 
strain, but independent of the strain itself, as stated by Newton's 
Viscosity Law. [8] At low shear stresses most liquids and gases obey 
Newton's Law: 



(2.2) T = 

where *T = applied shear stress 
0 = shear strain rate 
Y\ = Newtonian viscosity 
coefficient 

The laws and equations above represent idealizations. If sufficiently 
precise measurements are made it is found that no real solid is truly 
Hookean and that no real liquid is exactly Newtonian in mechanical re- 
sponse. [9] To represent such deviations an equation can be written in 
the form: 



Y = BT t'*" 



For the ideal elastic solid, /J'L equals zero and the equation becomes: 



(2.4) 



Y = BJ 
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For 



where the constant B is the reciprocal of the shear modulus, 
a perfectly viscous liquid, the exponent /YL equals one, and the equation 
becomes : 

r = art 

or 

(2.5) T — Bt 

where the constant B is the reciprocal of the Newtonian Viscosity, 

For a real material the exponent /Yl can vary between zero and one, and 
is a measure of the relative importance of the elastic and viscous contri- 
butions to the mechanical behavior. For values of the exponent /7U be- 
tween zero and one, equation (2.3) indicates a time behavior in which 
stress depends on both strain and strain rate. A material such as this 
which shows characteristics intermediate between those of an ideal solid 
or ideal liquid is called viscoelastic. 

In some situations neither the strain nor the strain rate is direct- 
ly proportional to stress. An example of this complicated relationship 
is the plastic deformation of metals. In the absence of this complexity, 
and when both strain and strain rate are linear with stress as indicated 
by equation (2.3) a material is said to display linear viscoelastic be- 
havior. [2] 

Linear viscoelasticity might be expected with materials whose struc- 
ture is not changing during the time it is subjected to an applied stress. 
[11] For such linear viscoelastic behavior, a law known as the Boltzmann 
Superposition Principle applies. This principle is illustrated in Figure 
1, which relates to a creep test where a material is subjected to a con- 
stant stress, and the creep strain changes with time. For linear visco- 
elastic behavior, this strain is directly proportional to stress: 
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< 2 . 6 ) e(t) = o 3(t) 

where T("t) > the creep compliance, is a function representing the time 

effects of the viscoelastic response. For a stress J 0 applied at zero 
time, the strain as a function of time is: 



(2.7) e(t) = <7 0 J(t) 

For some higher stress (X 0 + applied at zero time, the strain would be 
increased : 



(2.8) £(t) = (cT 0 -KT,) * J(0 

Now if (Jo i s applied at time "t = O , and stress CT, is added at a 
later time t - 1, the strain becomes: 

( 2 . 9 ) e(t) = o;- J(t) + <r, -J(t-t,) 

This is illustrated in Figure 1 where linear superposition of the strains 
at their elapsed times is indicated. The stress history may be general- 
ized by: 

tl - 1 

( 2 .10) e(t) = S <r r 3(t-U) 

tl--Oo 

For a continuous stress history starting at zero time, the Boltzmann 
Superposition Principle can be expressed in integral form: 

0- t 

3(t-e)or(e) 

=o 




( 2 . 11 ) 



e(t) = j 
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Stress 
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The Boltzmann Superposition Principle Applied to Creep (from Ref. 15) 



where Q is current time during the loading history, where cr(e) is 
written for > and is con stant in the integration. The creep 

strain in equation (2.11) is therefore represented by a convolution in- 
tegral where (j(Q ) is the "forcing function," and 3 (*t is the "re- 
sponse function. " 

In a stress relaxation test, a tensile strain is rapidly applied and 
the stress changes with time. For linear viscoelastic behavior. 



(2.12) cr(t) = e-E r (t) 

Here the viscoelastic time behavior is given by , called the re- 

laxation modulus . Equations analogous to (2.7) through (2.10) can be 
written for successively applied strains, so that for a continuous strain 
history, 



9=t 



( 2 .i 3 ) c r(t) = f Er(t-e) e(e~) 

e - o 



de 



It has been shown by Volterra that equations (2.11) and (2.13) are 
consequences of each other. [11, 12, 13] These equations are also a 
statement of Boltzmann's Superposition Principle and a definition of 
linear viscoelastic behavior. A thorough treatment of the mathematical 
theory and implications of linear viscoelasticity can be found in refer- 
ence [14]. 
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3. Linear Viscoelastic Models . 

Mechanical models are not necessary to describe viscoelastic phenom- 
ena, but they have been useful in describing a combination of Hookean and 
Newtonian characteristics of materials* The elastic mode of stress re- 
sponse may be represented by a spring (Figure 2), presumed to behave in 
an ideal Hookean manner, i.e., strain proportional to applied stress and 
independent of strain rate. The flow mechanism is indicated by a dashpot, 
(Figure 3) assumed to be filled with a Newtonian fluid so that its de- 
formation is proportional to the time and to applied stress, and inversely 
proportional to the viscosity of the fluid, . The behavior of a 

material displaying flow superimposed upon elasticity can be represented 
by a combination of a dashpot and a spring. The two basic models are the 
Maxwell Element (Figure 4) and the Voigt Element (Figure 5). The response 
of these units to an applied stress is illustrated in Figure 6. In order 
to represent a real material an infinite number of these elements may be 
required. Generalized Voigt and Maxwell models are shown in Figures 7 and 
8. Poincare [16] indicates that a phenomenon representable by any one 
model is also representable by an infinite number of other models. If 
one dashpot in the generalized Maxwell model has infinite viscosity 
( - 00 ) and if one dashpot in the generalized Voigt model has zero 

viscosity ( 7 \ = O ) then the models become equivalent in their representa- 
tion of viscoelastic behavior. 

The creep behavior of most materials can be considered as the super- 
position of three stress responses. (Figure 9) The instantaneous elastic 
strain can be represented by a Hookean spring, and the non-recoverable 
viscous flow by a Newtonian dashpot. The Voigt Element is particularly 
useful in describing a delayed elastic response. The viscous element 
coupled in parallel with the elastic element represents a damping resistance 
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Figure 2 
Hookean Solid 
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F igure 3 

Newtonian Liquid 
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Figure 4 
Maxwell Element 




Figure 5 
Voigt Element 
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Figure 6 

Response of Viscoelastic Models to Stress 

(a) Hookean Solid 

(b) Newtonian Liquid 

(c) Maxwell Element 

(d) Voigt Element 
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The Generalized Maxwell Model 





t 





Figure 9 

Creep Behavior of a Viscoelastic Material 
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Figure 10 

The Four-Element Model For Creep Behavior 
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to the establishment of elastic equilibrium. [10] A corresponding mech- 



from their minimum energy unstressed positions, which involves diffusion 
of chain segments and takes considerable time. From the equations re- 
presenting the components of the Voigt element, 



The solution to this governing differential equation gives the strain 
time relation: 



The creep behavior of linear viscoelastic materials illustrated in 
Figure 9 may then be described by the four-element model of Figure 10. 
The resulting equation becomes: 



The four-element model has been used successfully by Gearhart and 
Kennedy [17] in their study of cellulose acetate butyrate plastics, and 
by Dienes [18] to correlate data on creep of polyethylene at elevated 
temperatures. Marin and Pao [19] have extended the four parameter equa- 
tion to account for non-linear stress behavior frequently encountered in 



anism in a polymeric solid might be the uncoiling of long polymer chains 





(3.3) £ 




where 




is the retardation time, the time required for 



the delayed strain to reach a value 




9 



or 63.2% of its final 



value. 
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the delayed strain and viscous flow components for a simple plastic. 
Their equation (in its simplest form) can be written: 



( 3 . 4 ) £ = 



A modified version of this equation was developed to characterize the 
creep behavior of the materials studied in this investigation. 
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4. 



Non-Linear Viscoelastic Behavior. 



In a previous study [6] involving polymeric materials in the same 
loading ranges used in this investigation, marked non-linear behavior was 
encountered in the elastic, delayed elastic, and flow regions of stress 
response. Such deviations from linear viscoelasticity are not uncommon. 

[20] In a study of creep in cellulose acetate, Haward [21] found that the 
delayed elastic strain component increased by a factor of six when the 
stress was increased 30%. This non-linear behavior occurred at a compara- 
tively low stress of 500 psi. Similar stress non-linearities have been 
observed for glass-like amorphous polymers [22], polymers in the transi- 
tion region just below T [2], and polymers in the rubbery state [23]. 

S 

In all cases where deviations from linearity occur, an increase in stress 
produces a more- than-proportionate increase in strain. 

Alfrey [10] demonstrated that for large values of strain, non-linear- 
ity can be expected. From the theory of elasticity the work required to 
distort a solid to a given strain, called the "free energy of deformation, " 
can be represented by a power series: 

(4.1) tf'(Y) - + + C 4^ 4 

The shear modulus G has a magnitude given by: 

d z 7 

(4.2) G - d yt. 

For small values of strain only the quadratic term is significant and 
the modulus is a constant, independent of strain . At larger strains, 

when higher order terms become important, the shear modulus becomes a func- 
tion of , and stress non-linearities result. A similar analysis 
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predicts non-Newtonian flow of viscoelastic materials at higher stresses. 

For glassy polymers which have high moduli and can support large 
stresses, non-linear behavior has been observed even when the strains are 
relatively small. [2] This suggests that internal structural changes are 
responsible for the deviations. For example, it is reasonable to suppose 
that the flow rate will be markedly different for polymer chains which 
are straight than for a coiled chain. Experimental measurements in this 
field are sparse, but some theoretical aspects have been discussed by 
Schwarz 1 in terms of transition-state kinetics. [24] 

Nonlinearity of viscous and elastic elements of a model create 
systems of nonlinear differential equations which are difficult to solve 
by ordinary methods. In this investigation an equation similar to that 
of Marin and Pao (Equation 3.4) was used to fit creep data for various 
polymers. By making some simplifying assumptions in the equation, and 
with the aid of a digital computer, satisfactory predictions of composite 
behavior were obtained. 
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5. 



The Experiment. 



The viscoelastic materials selected for the investigation were: 
high density linear polyethylene, isotactic polypropylene, nylon 6/10, 
and polycarbonate. All were obtained from commercial sources in the form 
of extruded sheet. General mechanical properties of these plastics deter- 
mined by tension testing, along with approximate values of other physical 
constants from the literature [15, 25] are tabulated in Appendix A. 

Creep test specimens were cut and machined to the dimensions shown in 
Figure 11. The thickness of the specimens was that of the sheet materi- 
al; approximately 0.060 inches for polycarbonate and nylon, and 0.192 inches 
for polypropylene and polyethylene. Specimens used for relaxation tests 
and tensile tests were made with the dimensions shown in Figure 12. 

All specimens were cut with the long axis parallel to the extrusion direc- 
tion. 

Creep tests were carried out on constant load creep machines in a 
laboratory with the temperature maintained at 23 + 1° Celsius and relative 
humidity controlled at 50 + 2%. Construction features of the machines 
and the method of gripping the specimen are discussed in Appendix G. 

Series of creep tests were run at different stress levels, ranging be- 
tween 10 and 80 percent of the nominal yield strength. Extension of the 
test specimens over a two-inch gauge length was measured with a Tripoli- 
tis type extensometer which could be read to + 0.0001 inches. After load- 
ing the specimens, readings were taken at elapsed times of 6, 12, 30, 60, 

90, 120 seconds, and then in increasing intervals from 1 to 10 minutes 
during the first hour, then every 30 minutes during the first day. The 
test was continued for 100 to 150 hours, with an average of 50 data 
points taken during the run. Similar test procedures were followed for 
the composites, with the component materials clamped together outside the 
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Figure 12 

Tensile Test Specimen 




Figure 13 
Clamped Composite 
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gauge length to insure uniform elongation between the clamps. (See Figure 
13 ) 

Uniaxial tension tests were run at the same temperature and humidity 
on an Instron testing machine. Specimens were pulled at a constant rate 
of deformation, and Tripolitis extensometers were used to measure the 
extension over a two-inch gauge length. Preliminary experiments measur- 
ing the strain at timed intervals indicated that the strain rate was es- 
sentially constant throughout the test. 

Stress relaxation tests were performed on the Instron tester in the 
following manner. Samples were strained rapidly at ten inches per minute 
to a predetermined load level, where a load limit microswitch automati- 
cally stopped crosshead motion. Graphical output of load relaxation 
versus time at constant strain was obtained from the recorder installed 
on the testing machine. Initial straining was essentially instantaneous; 
less than one second was required to reach the desired stress. Dura- 
tion of the relaxation tests was approximately thirty hours. 

A portion of the creep data for polypropylene and polyethylene was 
obtained in a previous study. Because this work [6] was done in the 
same laboratory utilizing identical test procedures and specimens cut 
from the same stock, no distinction is made in the data tabulation in 
Appendix A. 
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6. Selection of an Equation to Characterize Creep Behavior. 

The four-element model (Figure 10) discussed in Section 3 describes 
the creep process in three stages; an instantaneous elastic strain, a 
delayed elastic response, and non-recoverable plastic flow. While there 
is considerable theoretical justification for this description of creep 
behavior [11, 27], this simple form of the model has two serious short- 
comings. First, it fails to account for the non-linear phenomena fre- 
quently encountered in tests made at stresses above 25% of the yield 
strength. Improved fits with actual data can be achieved by modifying 
the basic equation (Equation (3.3)) to account for these stress anomalies. 
Such a modification is shown in equation (3.4), and this method was uti- 
lized for data fitting in this investigation. 



The second limitation imposed by the four-element model is more funda- 
mental. The presence of only one Voigt Element implies that the delayed 
elastic response of a polymer is a single molecular mechanism with one 
characteristic retardation time. As pointed out by Alfrey [15], the coils 
and linkages of a polymer chain have a multitude of sizes and spatial ar- 
rangements. One would expect local kinks to be straightened out more 
rapidly than long range coils, and that elastic equilibrium would be 
achieved sooner if the chain segment had a favorable orientation with^the 
stress axis. The delayed elastic response of a viscoelastic material 
would then consist of several mechanisms with different retardation times. 
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The 55 memory" effects observed in high polymers are frequently cited as 
proof of the existence of more than one retardation time. Failure of a 
single Voigt Element to represent accurately the delayed elastic response 
in the materials used in this study is discussed in Appendix F. 

Figure 14 illustrates schematically the shape of a creep curve for 
a viscoelastic substance compared to a set of simple exponential curves 
of the form € - [*C ( \ — JL ^ ) , and shows the failure of a single 

time constant concept to represent actual data. Figure 15 shows the re- 
tardation spectra corresponding to the simple exponentials. In this figure, 
the term \d represents the rate constant (= Vr ) and the rate 

constant distribution function. The creep curve of Figure 14 can be ap- 
proximated very closely by a finite number of simple exponentials appropri- 
ately weighted. The model representing the delayed elastic response would 
then be the generalized Voigt Model of Figure 7. Such a representation 
consisting of exponentials CD , (^) , and (2) is given in Figure 16. 

A generalized equation for /H^ discreet relaxation times, with the in- 
stantaneous strain (A^) and the flow strain (A^ x t) included would be: 



( 6 . 1 ) 



-bit 



( ✓a. —Q* L- \ 

\ — 2 CL i fL* J 



where CLl = ^(b i) , the discontinuous distribution function for the 
exponential, and 2,^4= 1. Although this "spike distribution" model 
was not used in this study to correlate creep data, G. F. Kinney, in a 
parallel investigation at NPGS, has used the model to fit the same data. 
In his work, rate constants were arbitrarily assigned, then appropriately 
weighted using a least squares curve fitting technique adapted to a digi- 
tal computer. Two to six retardation times were used, and excellent fits 
were obtained, even with the simple two-spike distribution. 
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Figure 15 

Retardation Spectra Corresponding to 
the Simple Exponential Curves of Fig. 14 
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Figure 16 

Representation of Creep Behavior By Spike 
Distribution of Retardation Times 
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Another method of representing delayed elastic behavior is by a 
continuous distribution of retardation times, ( b ) . This concept is 
attractive from a theoretical standpoint if one considers delayed elasti- 
city as consisting of many molecular mechanisms whose rate constants form 
a quasi -continuous spectrum. Alfrey [10] develops a molecular theory of 
delayed response and arrives at a theoretical distribution of retardation 
times for a monodisperse amorphous linear polymer. This distribution, 
whose upper and lower limits can be determined from the theory, is shown 
in Figure 17. (Ihe notation here is left in the form given by Alfrey) 

For a polymer consisting of two different molecular weights, Alfrey has 
deduced the retardation spectrum shown in Figure 18. Tobolsky, Andrews, 
and Dunne 1 [28] used a box distribution (Figure 19), to correlate stress 
relaxation data. The limits of such a box distribution may be obtained 
graphically as shown in Figure 19 from a data plot of delayed elastic 
strain versus logarithm of time. 

The calculation of the retardation spectrum for an actual material, 
where the time constants extend over many decades of time, can be an in- 
volved process. A simple method [15] of approximating the spectrum is 
given by: 

(6.2) L(T) = 



Here,[ 3 (4.) — ] corresponds to the delayed strain response, 

^ 1— is the distribution of retardation times which 

can be converted to the rate constant distribution <}(b) . Figure 20 
demonstrates this graphically by plotting the slope of the delayed strain 
vs. logtime curve versus logtime. It is apparent that if Figure 20 does 
represent the response of the four-element model of Figure 10, the approxi- 
mation gives a cun/e with a peak at the most probable ^instead of a single 
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Figure 17 

Retardation Spectrum for a Monodisperse 
Amorphous Linear Polymer (Ref, [10] ) 
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Spectrum for a Polymer Consisting of 
Two Different Molecular Weight Fractions 
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Figure 19 

Box Distribution Of Retardation Times 
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Figure 20 

Method Of Approximating Retardation Spectra 
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spike. For real materials, however, the approximation is better. Refer- 
ence [27] reviews more accurate methods for estimating the retardation 
spectra. 

Analysis of the box distribution function is given in Appendix F. 

It is shown that this distribution is unsuitable for characterizing the 
creep behavior of the materials tested, as are other monomodal distribu- 
tion functions. 

In the absence of a mathematically tractable function to represent 
the retardation spectrum of the materials tested, the continuous spectrum 
concept was discarded. Instead, a simple empirical equation was developed 
by assuming that the rate constant associated with the delayed elastic 
response is a function of time. One might expect that a polymer chain 
segment should exhibit a changing retardation time as the uncoiling pro- 
cess continues and the shape of the molecule changes. Thus, an overall 
characteristic time constant might be considered, with a simple 

power dependency on time : 



Substitution of this expression into the equation for the four element 
model (Equation (6.4)) leads to an equation of the form of equation (6.5) 



(6.3) 
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Equation (6,5) was used to fit creep data in this investigation. The 
method of obtaining the creep parameters A^, A A^, C, and Q are dis 
cussed in the next section. 
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7. 



Analysis of the Data. 



Creep data were taken for all specimens in the form of elongation 
versus time. These data were punched on cards and a print-out obtained 
from an IBM 360 Digital Computer. The printout served as a check to pick 
up gross errors due to equipment malfunction during the testing period or 
errors made in transcription of the data. 

The first step in fitting the creep data was to obtain the time- 
independent part of the elastic strain. This instantaneous elastic 
strain, , occurring when the load is applied, is represented by the 

parameter A| in equation (7.1): 

ct Q x 

c 7 * 1 ) £ = + A^t +• ) 

To evaluate the constant /K\ , a direct method was employed. The first 
four data points were used to perform an extrapolation back to zero time. 
This extrapolation to determine the parameter is illustrated in 

Figure 21. It is clear from this figure that this evaluation of At 
is at best a rough estimate due to the nonlinear nature of the extrapola- 
tion. It will be shown in Appendix D that the parameter Aj does in- 
clude a portion of the delayed elastic strain, but for the time frame of 
this experiment, the strain At occurs essentially instantaneously. 

To account for the non-linear viscoelastic behavior encountered, the 
parameters A| > A^ , and A 3 were correlated with stress so that 
the creep behavior could be predicted for any stress level. The general- 
ized form of equation (7.1) becomes: 

<7 ' 2> e = Kct n + &cr R t + Po- M (\-^ ctS> ) 
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Cr<z<z.p Strain , in. /in x lO 3 




Figure 21 

Extrapolation of Creep Data to Determine Instantaneous Strain 
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The magnitude of parameters A?_ and A 3 depend on the value chosen for 
A( • Prior to evaluating A 2. and A 3 , the extrapolated values 
of zero time strain A| were correlated with the first term of equation 
(7.2): 

(7.3> a, = Kcr* 

( 7 . 4 ) A, = JUK +■ mVct 

Values of J/y c A\\ for different stress levels were read into a digital 
computer, and equation (7,4) was used to perform a least squares fit to 
determine the stress coefficient and the stress exponent H 

Figure 22 shows the plot of equation (7.4) and the fitted curve. The 
data displayed in this figure are that obtained from the extrapolations 
of Figure 21, made for polycarbonate. 

An independent method of obtaining the zero-time strain was attempted 
by running a series of tensile tests at various strain rates, and plotting 
the resulting data as sketched in Figure 23. From these stress-strain 
curves, a cross-plot of strain versus reciprocal strain rate may be con- 
structed as illustrated in Figure 24. The zero-time strain can then be 
obtained by extrapolating the data to strain rate £ = , which 

corresponds to instantaneous loading. This technique was attempted using 
a digital computer to perform a parabolic extrapolation, but the exponential 
character of the data plot invariably resulted in high estimates of the 
zero-time strain, £o , as indicated in Figure 24. In addition, the 
maximum deformation rate available on the Instron tester used for the 
tensile tests was £ = 10 inches /minute . Consequently, there was a lack 
of data in the region of greatest interest, i.e., where J/^ 0. This 
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Parameter Ai inchcs/inch 




Figure 22 

Correlation of Parameter A-^ with Stress 
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Figure 23 

Stress-Strain Curves for Various Strain Rates 




Figure 24 

Zero Time Strain Extrapolations from Stress-Strain Data 
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method may be expected to produce excellent results if high-speed testing 
apparatus were available [29], but in view of the limitations discussed 
above it was discarded in favor of the creep data extrapolation technique. 

After the parameter was found by extrapolation and ^smoothed" by 
the stress correlation (Equation (7.3)) the creep data, along with the 
calculated parameter A., was placed in the digital computer via program 

^ i 

FITDATA (See Appendix B) . The viscous flow parameter A ^ was determined 
by a linear extrapolation using data points at elapsed times greater than 
70 hours. It was assumed that elastic equilibrium was achieved by this 
time. With parameter A 1 specified and parameter A^ computed, the maximum 
magnitude of the delayed elastic strain, A,, is determined. This pro- 
cedure used by program FITDATA is outlined in the sketch of Figure 25. 

The time constant parameters were computed by FITDATA by evaluating the 
delayed strain fraction, , at various times. 

(7,5) <5 = ^/A-i = I " 

from which . ^ 

= ct ft 

(7-6) = Zn-C ■*" Q Zn/t 

A plot of the left-hand side of equation (7.6) versus was given 

a least squares curve fit, and the constants C and Q determined. The 
method is sketched in Figure 26, and graphical computer output of actual 
data is displayed in Appendix B. 

The creep parameters A^ and A^ were correlated with stress using 
equations (7.7) and (7.8): 

(7.7) A t = &cr R 
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Figure 25 

Determination of Parameters A 2 and A 3 



ln[-b(l-£)j 



t= I Hr., lnt*0 



Figure 26 

Determination of Parameters C and Q 




Int 
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(7.8) 



A, = Po- M 

The method used for the A correlation was employed to determine the stress 
exponents and coefficients in equations (7.7) and (7.8). These determina- 
tions for polycarbonate are illustrated in Figures 27 and 28. 

A question arises as to the validity of inclusion of the viscous 
flow term A^/t, in equation (7.1). If the term is not included and 
the first and second derivatives of strain with respect to log time are 
taken, at long times the second derivative is always negative and becomes ’* 
vanishingly small. At long times the first derivative, corresponding to 
the slope of the 6. vs. curve of Figure 29, approaches zero. With 

no viscous flow, elastic equilibrium is eventually reached and extension 
ceases. If the flow term t is included, the second derivative 

changes sign and becomes increasingly positive at long times. This behav- 
ior is observed as an inflection and an upswing in the plot of £ vs. , 

Figure 30. If viscous flow is a factor in creep behavior, therefore, a 
data plot of £ vs. "t should show the inflection and upturn indi- 
cated by Figure 30. In all the materials investigated this behavior due 
to viscous flow was observed, as illustrated by the example of Figure 31, 
justifying inclusion of the term A 2 . * "t in equation (7.1). 

Values of the creep parameters determined for the materials tested 
are summarized in Table 2 of Appendix A. The constants C and Q are ap- 
parently independent of stress and are of the same magnitude for all 
materials tested. Standard deviations in strain were acceptable for all 
runs. The errors in computed strain were less than 5% in the region where 
delayed elastic behavior is significant, except at elapsed times under 
two minutes. Sample computer output from FITDATA listing point by point 
comparison of observed and computed strain is given in Appendix B. For 
typical values of the constants C “ 0.8 , and Gt - 0.85 , the delayed 
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Parameter Ai , in/in/hr. x 10 




Figure 27 

Correlation of Parameter A£ with Stress 



50 




Parameter A3 , in /in x 10 
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Figure 28 

Correlation of Parameter Ag with Stress 
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Figure 30 

Strain versus Ln Time with Viscous Flow 



52 




53 



time , 



strain is only 97% exhausted at an elapsed time of 70 hours. For the 
extrapolation to determine parameter elastic equilibrium is assumed. 

Consequently, in the resulting "forced fit, 1 ' after the other parameters 
are determined, the computed strain is approximately one to two percent 
low at long times. 

From the tests at higher stresses, it is apparent that elastic equi- 
librium is not achieved during the testing period. This produces a high 
estimate of the parameter A^, which in turn results in a low estimate for 
parameter A^* At low stresses the flow rates are very small and the er- 
ror introduced is negligible. However, power dependence of A^ on stress 
creates noticeable deviations at high stress levels. This deviation is 
clear in Figures 27 and 28 at stress levels of CT = 4500 psi and Cf = 
5000 psi. In performing the stress correlations, data at lower stresses 
was weighted more heavily, which corrects the erroneous values of A^ and 
A^. (See figures 27 and 28). 

The fitted empirical equation (7.2) provides excellent predictions 
of creep behavior on individual materials loaded at constant stress. 

When a simple composite sandwich of two materials is tested in creep how- 
ever, the stress in each component changes with time due to the differ- 
ent viscoelastic responses of each material. If the time responses were 
linear in stress, accurate theoretical predictions of composite behavior 
could be made. Reference [15] outlines procedures for calculating stress 
as a function of strain and time when the creep compliance j(t') is 
linear in stress. The Boltzmann Superposition Principle (Equation (2.11)) 
could then be applied to determine the creep behavior of a composite. 

From Table 2 in Appendix A, the stress exponents are greater than one, so 
the viscoelastic behavior of the materials tested is non-linear and the 
creep compliance is not independent of stress. Resulting stress- 

strain-time relationships become extremely complex. If a non-linear 
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viscoelastic "equation of state" could be developed, it might be too compli 
cated for practical use, even with digital computer techniques. 

In order to predict the behavior of composites studied in this investi 
gation, a simplifying assumption was made. The strain response of a 
material was assumed to be independent of its stress-time history. Equa- 
tion (7.2) then functions as an equation of state; if a stress and time 
are specified, then the strain is determined. From the physical restric- 
tions placed on the composite tested in creep, 

(7.9) P*u -V- ?b - P 

(7.10) €«. = €-b ~ € 

where P represents applied load and the subscripts a and b are used to 
identify quantities pertaining to component A or component B of the com- 
posite. Letting a and b denote the relative fraction of component A 
and component B respectively, 

(7.11) GL + b = t 

(7.12) GL01 + b0* b = cr 

Unsubscripted symbols above pertain to the composite structure as a whole. 
From equation (7.2), 

xQb 

(7.13) e b = + Pbcr b b (\-/ ) 

< 714 > + B*cr^t + ?^(\ 
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Equation (7.12) may be solved for Gcu in terms of known composite stress 
(X , known volume fractions a and b, and the unknown stress G± . 

With this result substituted for Go- in equation (7.14), equations (7.13) 
and (7.14) may be equated according to equation (7.10). For a specified 
time, the corresponding value of Gb may be found. A series of values 
of Gb and time can then be substituted in equation (7.13) to determine 
strain €~b according to the equation of state assumption. By equation 
(7.10), strain is the same as the strain in the composite. 

A digital computer program, program COMCUR, was used to predict compos 
ite behavior in this manner. Sample output from COMCUR consisting of point 
by-point comparisons of predicted and observed strain is found in Appendix 
C, along with graphical displays of this output. A summary of results of 
this prediction method are tabulated in Table 4 of Appendix A. When the 
viscoelastic responses are similar, the load distribution in the components 
of the composite does not undergo large changes during the test and pre- 
dictions are excellent. Such is the case for the polyethylene/polypro- 
pylene (PE/PP) composites. It should be noted that if the stress in the 
components were constant, the method would be theoretically correct, since 
equations (7.13) and (7.14) are valid for a constant stress test. Large 
errors in predictions occur with component materials of nylon and poly- 
carbonate. As indicated by the magnitude of the creep parameters given 
in Table 2, Appendix A, the viscoelastic response of these materials is 
quite different. Graphical output from COMCUR shows that the load dis- 
tribution for the NY/PC composites changes rapidly during the first few 
minutes of the test. 

The empirical creep equation (Equation 7.2) was used to predict stress 
strain behavior and stress relaxation behavior for single materials. This 
analysis was done with a digital computer using the equation of state 
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assumption. An explanation of program RELAX, used to predict stress re- 
laxation behavior is given in Appendix D. Output from program SIGEPS, an 
estimate of stress-strain curves, is given in Appendix E. In both cases, 
computed responses were considerably different from data obtained by experi- 
ment . 
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8. Conclusions. 



The results of computer programs RELAX and SIGEPS demonstrate the 
complex nature of viscoelastic behavior. In tensile tests and relaxation 
experiments the stress changes with time and the constant stress equation 
gives poor predictions of responses. While the equation gives excellent 
estimates of creep behavior, the stress-time history cannot be ignored 
when predicting viscoelastic behavior in other modes of testing. 

The combined mathematical and experimental approach to understand- 
ing composite behavior is a powerful tool. An analysis of simple compos- 
ite structures provides guidelines for predicting behavior of more compli- 
cated systems where factors such as component distribution and interfacial 
shear stresses become important. The non-linear behavior of many visco- 
elastic materials, even at moderate stresses, leads to involved mathemati- 
cal relationships. A digital computer is necessary for analysis of these 
complex systems. 

Some opportunities for further investigation are: 

(1) Use of X-ray diffraction techniques to follow structural 
changes during testing in order to determine the relation 
of internal structure change to observed non-linear visco- 
elastic effects. 

(2) Fabrication and testing of more realistic composite materi- 
als to determine if the empirical equation developed here is 
general enough to characterize the behavior of these materi- 
als in creep. 

(3) Use of other modes of testing, such as stress relaxation and 
dynamic mechanical tests to assist in development of non- 
linear viscoelastic theory. 



58 



BIBLIOGRAPHY 



1. Dietz, A. G. H. , Composite Materials , American Society for Testing 

and Materials, Philadelphia, Pa., 1965. 

2. Ferry, J. D. , Viscoelastic Properties of Polymers , John Wiley and 

Sons, Inc., New York, 1961. 

3. Special Report, Materials in Design Engineering , Vol. 16, No. 4, 

1965. 

4. Holliday, Composite Materials , Elsevier Publishing Co., 1966. 

5. Schapery, R. A., "Stress Analysis of Viscoelastic Composite Materials", 

Journal of Composite Materials , Vol. 1, No. 3, pp. 228-266, 1967. 

6. Howard, M. A., Thesis: "Constant Temperature Creep Studies of Some 

Simple Polymeric Composites ", Naval Postgraduate School, Monterey, 
Calif., June 1967. 

7. 1965 Book of ASTM Standards - Part 27, Plastics, General Methods for 

Testing, American Society for Testing and Materials, Philadelphia, 
Pa., 1965. 

8. Andrade, E. N. da C. , Viscosity and Plasticity , Chem. Publishing Co., 

Inc. , New York, 1951. 

9. Zener, C., Elasticity and Anelasticity of Metals , University of 

Chicago Press, Chicago, 111., 1948. 

10. Alfrey, T. Jr., Mechanical Behavior of High Polymers , Interscience 

Publishers, Inc., New York, 1948. 

11. Tobolsky, A. V., Properties and Structure of Polymers , John Wiley 

and Sons, Inc., New York, I960. 

12. Volterra, V., Theory of Functions , Dover Publishers, Inc., New York, 

N. Y., 1959. 

13. Leaderman, H. , Elastic and Creep Properties of Filamentous Materials 

and Other High Polymers , The Textile Foundation, Inc., Washington, 

D. C., 1943. 

14. Gross, B., Mathematical Structure of the Theories of Viscoelasticity , 

Hermann et Cie, Paris, 1953. 

15. Nielsen, L. E., Mechanical Properties of Polymers , Reinhold Pub. 

Corp. , 1967 . 

16. Poincare', H. , The Foundation of Science , The Science Press, New York, 

1929, p 181. 



59 



17. 



Gearhard, W. M. and Kennedy, W. D., Ind. Eng. Chem . , Vol 41, p 695 
ff., 1949, 

18. Dienes, G. J. , Journal of Colloid Science , Vol 2, p 131, 1947. 

19. Pao, Y. and Marin, J. , Journal of Applied Mechanics , Vol. 20, 

No. 2, pp. 245-252, 1953. 

20. Reiner, M. , Rheology , Vol. 1, Ch 2, Eirich, F. R. , Ed., Academic 

Press, Inc., New York, N. Y. , 1956. 

21. Haward, R. N. , Transactions, Faraday Society , Vol. 38, p 394-403, 

1942. 

22. Lethersich, W. , British Journal of Applied Physics , Vol. 1, p 294, 

1950. 

23. Leaderman, H., Journal of Polymer Science , Vol. 14, p 261, 1955. 

24. Schwarzl, F., Kolloid-Zeitschrif t , Vol. 165, p 88, 1959. 

25. Kinney, G. F. , Engineering Properties and Applications of Plastics , 

John Wiley and Sons, Inc., New York, 1957. 

26. Keith, H. H. , "Instrumentation for Mechanical Testing *', (Notes from 

faculty seminar, U. S. Naval Academy, Annapolis, Md., May, 1965.) 

27. Leaderman,' H. , Rheology , Vol. 2, Ch 1, Eirich, F. R. , Ed., Academic 

Press, Inc., New York, N.Y., 1956. 

28. Tobolsky, A. V., Dunell, B. A., and Andrews, R. D. , Textile Research 

Journal , Vol. 21, 1951. 

29. Deitz, A. G. H. , and Eirich, F. R. , High Speed Testing, Vol. I , 

Interscience Publishers, Inc., New York, 1958. 

30. Pao, Y., and Marin, J. , Proc. American Society for Testing and 

Materials, Philadelphia, Pa., Vol. 2, 1952. 



60 



APPENDIX A 



TABULATION OF RESULTS 



Table 1 

Properties of Tested Polymers 



Material 


^Tensile 

Strength 

psi 


*Tensile 
Modulus „ 
psi x 10 


*Percent 

Elong. 


Melt. 

Point 

°C 


Glass 

Transition 

T , °C. 

g 


POLYETHYLENE 


4000 


220 


380 


137 


-120 


POLYPROPYLENE 


4400 


190 


210 


176 


-10 


POLYCARBONATE 


8800 


300 


60 


267 


150 


NYLON 


5900 


120 


185 


227 


40 



*At strain rate = 0.012 in/in/min 



/<5 
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Creep Parameters at Various Stress Levels - Single Specimens 
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Correlated Stress Exponents and Coefficients 
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APPENDIX B 



COMPUTER PROGRAM FOR DETERMINING CREEP PARAMETERS 
A Fortran IV computer program, compatible with the IBM 360 system, 
was used to determine the creep parameters A^, A^, C, and Q in the 
equation: ^ 

€ — At + A^'t + — ) 

The program, called FITDATA, is entered with the predetermined parameter 
A ^ and the def ormation-time data obtained from the creep tests. Computer 
output consists of the calculated parameters, tabulated values of observed 
and computed strain, and intermediate data used in the parametric determina- 
tion. An off-line plotting routine, called SUBROUTINE DRAW, is used to 
display the output (experimental data and fitted curve) in graphical form. 
This subroutine was also used to plot the data for determining parameters 
C and Q in the equation above. Least squares fitting was done by the 
library subroutine, SUBROUTINE LSQPOL available at the NPGS Computer 
Facility. Program FITDATA, along with typical tabulated and graphical 
output, is given on the following pages. The program, discussed briefly 
in Section 7, is annotated with comment cards. 



67 



. COMPUTER PROGRAM “FITDATA" 



f UK TRAN IV G LFVEL I, HUD 1 



D AT F = 6S119 



OOCl 
0 0 C 2 
GC03 



OUCfc 
UCC 7 
uCu* 

uccs 

CC1C 

m\ 

C C 1 3 
OC 1** 
CC15 
uClfc 
SCI 7 
CC1P 
OC IS 

3:‘i 



^c‘i 



CC3e 

„Cj7 

uC;6 



IMPLICIT REAL»ti (A-H). REAL *6 I O-Z ) 

Rfc AL*H ITITLFI12I, KTITLEU2I. LARA/BH / 

DIMENSION kP(lOJ) ,1(100) , TZ ( 100) ,EPZ( 100 ) .Y! 1001 , DFL Y ( 1 

100) . B( 5) , $fl!5) , T0( 5) .ST (5) .CQl 51 »SC(5) , AQ(30. 30) ,F0( 100) ,TLN( 100) 
2.E0Z! 100UTLN7 (100) ,TGG(200) .YOU 100) ,FRR (100) ,FRRF( 100) , FFD( 1001 . 
3EP0( 100) »fn< 100) .A( 30. JOT.TTI 100 ).X 1100 ) , F 2 < 100) .Cl *5) , ED/ / ( IOC) 

*f AL LABEL ( *>0) ,LABH.TG(20U) , YG< 200) , YGm<2oO) » TS ( 100 ) , E ( 100) ,TA( IOC 
1 ), fcM( 100 l.Dn-rT.OIFEP.OIST.TSKAL.F PSKAL.fi, ML CH/<*H /, 

2PMKP/AH P/ .PM ARK ( 100) , T Y{ 10GJ , TINY! 100) ,TYS < 100) .EOS! 1001 ,TLO , 

3FLAG.CZ .ZSCAL.YSCAL ,TU< 100) , YU( 100) ,EDLO, W!100) 

REAL AA1,«A2»Aa3»CC(;,000 

FOR EACH COMPUTER RUN THE FOLLOWING INFORMATION IS R F AH IN 
NOAT = NUMBER OE SFTS OF DATA 

IX = NUMBER OF SETS OF OATA TO BF ORAWN ON FACH GRAPH. 
JC = CUUNTING SYMBOL TO CHANGE THE POINT PLOTTING 
SYMBOL FOR EACH OATA SET. 

LABEL = RUN NUMBERS USED TO LABEL EACH DRAWN CREEP 
CURVE. 

ITITLE = TITLE OF GRAPH OF CREEP CURVES 

KTITLF TITLE E OR GRAPH SHOWING C AND 0 DETERMINATION. 

K E AC 1 6 , NDAT.IX 

16 FUKMATI2I5) 

READ 17, (LABEL! I) , I = l.NDAT) 

17 FORMAT ( 20A4 ) 

READ 1 H , I ITITLtI I ) , I = 1,12) 
lfl F0RMATI6AB) 

NSKIP = 0 

READ 502, IKTITLEI I), I = 1,12) 

532 F0KMATI6AB) 

SC 3 A I X = IX 

FNCAT = NOAT 
PAGE - FNDAT/A I X 
NPAGE * PAGE 
NOO = NPAGE*IX 
kC - C 

OC *.1 K = l.NPAGE 

A NEW GRAPH SHEET IS STARTEO EACH TIME THE INOFX OF THIS 
00 LOOP IS CHANGED. 



*.iU2 




JC = 1 


cC2i 


*0 


OC 42 L = 1.IX 


uc;<» 




JC = JC *1 


DC 2 5 




KC = KC + l 


c C26 




PRINT 10 


0C2 7 


10 


FORMAT l 1 H 1 ) 



THE FOLLOWING OATA IS REAn IN: 

NPTS = NUMBER CF POINTS IN THE DATA SET. 

GL = TEST GUAGE LENGTH 
A 1 = PREDETERMINED ZERO TIME STRAIN 

TMAX = MAXIMUM TIME IN HOURS TO WHICH THE FITTED CURVE 
IS TO BE OPAwN. 

EPSKAL = STRAIN SCALE ON CREEP CUPVE GRAPH (IN/IN/INCH) 
TSKAL = TIME SCALE ON CREEP CURVE GRAPH (HOUPS/INCH) 
CONST = CONSTANT USEO TO CONVERT INPUT CREEP DATA 
TO STRAIN 

FUN = RUN OESIGNATOW USED TD IDENTIFY THE DATA 
EP = ARRAY OF CREEP DEFORMATION DATA; USUALLY READ Iw 
AS OEFORMAT ION IN MILS 

TT * ARRAY OF T 1 MF S CORRESPONDING TO TH«= APPAY "FP« 

REaC 23, NPTS, GL.A1 ,TMAX, EPSKAL , TSKAL, C INST, RUN 
2 ) FORMAT! 15, 6F10.C,5X,A5) 

READ 25, (EPII ) , I = 1,NPTS ) 

25 FORMAT! UF*». 1) 

RE AD 26, ( TT ( I ) , I = 1 ,NPTS) 

26 FORMAT! 7F1J. 5) 

CONVERT THE OEF OR MA T I ON DATA TO STRAIN IN TNC HE S / I^CH , 

MAKE THE INITIAL VALUE OF STRAIN FOUAL TO THE 7 FRE T I me 
STRAIN AND THE INITIAL TIME FOUAL ZERO. NUMBFR OF DATA 
POINTS IS INCRFASFO BY ONE, BECOMES "NPTSl”. 

e F r I 1 ) = A 1 

Tull) = C • 

JO 4 3 I = l.NPTS 
CPCU ♦ 1) = EF(II/(CCNST*GLI 
■** n iuii = tti i ) 

NPTSl = NPTS ♦ l 
OC 30 0 I - l, NPTSl 
E P ( I ) = EP ( I ) 

30R TT! I ) = Till) 



THc NEXT SIXTEEN CAPPS SELECT ELAPSED TImfs GPFATE<? jhaf 
7C X OF THF TOTAL TESTING TImr. IF LFSS THAN FOUR POINTS 
A Kr OBT A [NEC, THC NEXT LONGEST TJmE IS SELECTED /NO THE 
PROCESS IS CONTINUED UNTIL FCIJK POINTS A R F OPT A I NFP . IE 
THE SHORTEST TIME IS LESS THAN 50* PE T <E totm FLAPSPD 

time, tie cata is deemeo unreliable for dftepmiming 

PARAMETER A2, and IS OISCARDFD. THE S c L EC T E C TIm.c FOR 
THt- EATrtAPi.LATICN AK F ST )R E D IN THE AcCAY "T7*» A \ D THF 
CuRfi E SFUNOI NG STRAINS ARC STOKF) IN ThE ARRAY ”F P/". 



T w = ». • 5 * T T ! NP T > I ) 
TE - C. 7*U INPTS1 ) 



lC4 5 




K A = C 




CCnt 


3C9 


KA = KA ♦ 1 




C w 4 / 




IFITT(KA) - Tf I 30<U31C, 


310 




310 


KA - KA-1 




CC4S 




KB - NPTSl - KA 




CC5C 




IF !KH - 4 ) 31C.31I.3ll 




C C 5 1 


311 


IF < TT IKA* 1 ) .IT. TM| GO T 


. 777 


oC32 




GC TC 312 




Oc = 3 


777 


PRINT 7 7e » RUN 




Jt;<i 


7 7d 


FORMAT!///, 1CX, »OATA FDR 
LEEP LXTRAPULATI JN* ,//) 


kUN 


« c 6 6 




GO TC 46 




c C 5o 


312 


03 313 I = 1,K3 




0057 




TZ < I ) = TT (KA+ I ) 




^065 


313 


FEZ l I ) * E P ! K A ♦ I ) 





UNRELIABLE FOR ST E ATY-S TATE <K 



THE EXTRAPOLATION IS DONE USING LIBRARY S'JHRCUTINF 
LSJPOL. THIS ROUTINE PERFORMS A FIRST DFGRPF PUl y- 
N; * I AL FIT TO THE SELECTED POINTS. f n^FF IC I t NT S OF THE 
FITTED STRAIGHT LINF ARE STCRFD IN THF ARRAY CALL FO "R", 
T h f SUBROUTINE PRINTS OUT INPUT DATA AND TEH- RFMILTS np 
THc FIT IN TABULAR FORM. THE VALUES OF THF CRPt-P 
PA^Amet-RS A2 AND A3 ARE THEN PFTceMiMrf) 



u JU 
-C 1 1 

W ZL2 



cALL LS‘. POL (KB, 1 ,0.0,0, SIGMA, TZ ,FP7, 
IF ( S I (,MA ,GT .0. O "> ) GO TO 777 
A 2 = B ! 2 ) 

A 3 = H ( 1) - A 1 



W.Y.DFl Y,B,SD,TO,f T,CO,Sf ,50) 



P AG C CC01 



68 



0C63 

0064 

0065 
0C66 
0067 



0C79 

CC80 

0081 

0C82 

0083 

0084 



CC85 

0Cd6 

0087 

0088 
CC85 
0C90 
CC91 
0C52 



0C 53 
CC54 



0C55 
C C 56 
CC57 
0058 
CC55 
C ICO 
C1C1 
0102 
C 103 



0 1 4 1 
0142 
C 1 43 
C 144 



0145 
C 146 
0147 
C 1 48 
0145 
C15C 
C 1 5 1 
C 1 52 



COMPUTE THE OELAYEO STRAIN FRACTION AS ARRAY "ED" AND 
NATURAL LOGARITHM OF TIME AS THE ARRAY "TLN" 



314 Eom =* 
TIN! 1 ) 
00 322 
322 TLN! I) 



♦ A2*TT( I > ) ) /A 3 



0068 




LQ * 10 


0C69 




IP * 0 


DC 70 
OC 71 


339 


If ( i Ip 


0072 


340 


EDZ I IP) 


0C7 3 




TLNZ! IP 


0074 




IF! 1 .-El 


0075 


341 


K PT S = 


0C76 




00 325 


CC77 




EOZZ! I) 


0078 


325 


EOZZ! I ) 



I = 1.NPTS1 
IEP( I J - ( A 1 
* 0.0 

I * 2.NPTS1 
= OLOG! TT < I I > 

COMPUTE THE DATA REQUIRED TO DETERMINE CONSTANTS C AND 0. 
DISCARD THE LAST TEN POINTS WHICH MAY GIVE ED VALUES 
GREATER THAN ONE» DUE TO THE "FORCED FIT" OF PARAMETER A2. 
ELIMINATE ANY VALUES IN THE INITIAL PART OF THE RUN WHFRF 
ED MAY BE LESS THAN ZERO BECAUSE OF AN OVERESTIMATE OF 

parameter ai. atore the results of thf calculation in thf 

ARRAYS "EOZZ" AND "TLNZ". 



• EDI I P ♦ LQ ) 

! * TLN! IP ♦ LD) 



ULOGI EOZZ I I ) ) 

USE SUBROUTINE LSDPOL TO PERFORM A LINEAR FIT TO THE DATA 
"EOZZ" AND "TLNZ". PRINT OUT THE DATA AND THE RESUITING 
FIT. 

CALL LSQPOLUPTS.l ,0,0,0, SIGMA, TLNZ , E OZ Z , W , Y , DEL Y , fl, SB , TO, S T , CO , 
ISC, AO) 

00 205 I = 1 , K PTS 

FED! I ) = 8(1) ♦ 8 ( 2 ) *TLNZ ! I ) 

205 ERRFII) = IIEOZZII) - FED I I ) ) * 1 00. ) /E OZ Z I I ) 

PRINT 71, RUN 

71 FORMAT! ///, 10X, »THE DATA USED TO F I NO "C" AND "0" FOR RUN »,A5, 

1' IS LISTED BELOW: *,///, 6X, ‘DATA PT. LN TIME LN ( -LN ( 1 . -ED ) ) 

2FITTE0 VALUES PCT .ERROR * , // ) 

DO 210 I = l.KPTS 
11= I ♦ LD 

210 PRINT 72, II, TLNZII),EDZZII),FED(I)i ERRFII) 

72 FORMAT! I12,F12.4,F13.4,F16.4,F13.2) 

IFISIGMA.GT.0.1) GO TO 891 

GO TO 893 

891 PRINT 892, SIGMA 

892 FORMAT! /// ,10X , 'SIGMA = *,F5.3,' UNRELIABLE FIT FOR FINDING C AND 
1 0',//) 

FROM THE COEFFICIENTS OF THE FITTED POLYNOMIAL, COMPUTE 
THE CREEP PARAMETERS C AND 0. 

893 Q = B ( 2 ) 

CC = OEXPIB! 1) ) 

THE NEXT NINE STATEMENTS ASSIGN 200 VALUES TO THE ARRAY OF 
TIMES "TGG" IN EQUAL INCREMENTS FROM ZERO TO TMAX. 

VALUES OF FITTED STRAIN ARE COMPUTED IN FUNCTION SUBROGRAM 
"EON" USING THE DETERMINED CREEP PAR AMFT FRS, AND STOREO IN 
THE ARRAY "YG". 

R = TMAX/200. 

TG ( 1 ) = 0.0 
TGG ( 1 ) = 0.0 
DO 44 I = 2,200 
TG! I } = TG! 1-1 ) ♦ R 

44 TGG! I ) = TG! I ) 

00 45 I - 1,200 

YGII) = E ON! A 1 . A2 * A 3 , CC , 0 , T GG ( I ) ) 

45 YGM! I ) ^ -YG! I I 

THE NEXT EIGHTEEN STATEMENTS REDUCE THE DATA POINTS SO 
THAT NO TWO DATA POINTS ARE CLOSER THAN 0.2 INCHES APART 
ON THE GRAPH SHEET 



01C4 




El 1) = EP 1 1 ) 


0 1C 5 




T A ( 1 ) = TTll) 


01C6 




JA= 1 


01C7 




IZ = 1 


01C8 




KR = 2 


0 1C5 


60 


JA = JA ♦ 1 


C11C 




IF! JA.GT.NPTS) GU TO 62 


cm 




TSCAL = TSKAL 


0112 




ESCAL = EPSKAL 


0113 




OIFFT = ITT! JA) - TT! IZ 


0114 




OIFEP = (EP! JA) - EP! IZ 


0115 




0 1 ST = S0RTI0IFFT**2 ♦ 


0116 




I F 1 01 ST - 0.2) 60,61,61 


0117 


61 


TA(KR) = TT(JA) 


one 




E ! KR ) = E P ( J A ) 


Cil5 




KR = KR ♦ 1 


0120 




IZ = JA 


C 1 2 1 




GO TO 60 



THE FOLLOWING STATEMENTS (DOWN TO STATEMENT 50) FURTHER 
REDUCE THE NUMBER OF DATA POINTS TO THIRTY OR LESS, AS 
REDUIRFD BY SUBROUTINE DRAW. 



0122 


62 


KR * KR -1 


0123 




IFIKR.LE.3C) GO TO 50 


0124 




IF1KR.LE.59) GO TO 63 


C 1 2 5 




KO = KR/20 


Cl<t 




KPLCT = KR/KO 


CU7 




DO 47 I = 1 »Kf LOT 


C 12 8 




E ( I ) = El I *K 0 ) 


0125 


47 


TAI I ) = T A ( I *K 0 ) 


0 1 3C 




KK = KPLCT 


0131 




GO TO 50 


C 1 3 2 


63 


KD = KP - 30 


C133 




KDP1 = KD ♦ 1 


0134 




DO 64 I = 1 , KO 


C 135 




T A ! I ) = TA( 1*2) 


0126 


64 


E ! I ) = E( 1*2) 


0137 




00 65 I = KDP1.3Q 


C 1 38 




TA ( I ) * TAI I ,KD) 


C 1 35 


65 


E! I ) = E! I+KD) 


014C 




KR = 30 



THE I NT EGE R "NSK I F" DETERMINES WHICH GRAPH PLOTS WILL BF 
EXECUTED 

NSK I P = 0: DRAW ONLY THE CREEP CURVES 

NSKIP = 1: DRAW ONLY THE CURVES FOR DETERMINING C ANO 0 



50 IFINSKlP.NE.il GO TO 818 
DO ail I = l.KPTS 
TYS! I ) = TLNZ! I > 

811 EOS! I ) = EDZZI I) 



THE FOLLOWING STATEMENTS COMPUTE GRAPH SCALES FOR PLOTTING 
THE C AND 0 DETERMINATION. THE NUMBER OF DATA POINTS 
ARE REDUCED TO THIRTY OR LESS. 



TLO = A8SITYS! 1 > ) 

I TLO = TLO * 10 • 

TLC = I TLO ♦ 1 
TLC = TLC/10. 

IF! TYS! U.GT.O.O) 
6DL0 = AH S ( E OS ! 1 ) ) 
IEDLO = EOLO* 1 0 . 
EDLO = IEDLO ♦ l 



TLO = -TLO 
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0153 

0154 

0155 

0156 

0157 

0158 

0159 

0160 

0161 

0162 

0163 

0164 

0165 

0166 

0167 

0168 

0 169 

0170 

0171 

P 

0174 

mi 

C 1 7 7 

0178 

Cl 79 

0100 

0181 

0102 

0 03 

0104 

0105 

0106 

0107 

0108 

0189 

0190 

0191 

C 192 

0193 

C 1 94 

0195 

0196 

0197 

C1S0 

0199 

0200 

0201 

02C2 

0203 

C2C4 

0205 

C2C6 

0207 

02C8 

0209 

02 1C 

0211 

0212 

0213 

C 2 14 

0215 

0216 

0217 

0218 

C 2 19 

C22C 

0221 

0222 

0223 

0224 

0225 

C226 

0227 

0220 

C229 

02 30 

C 2 3 1 

0232 

0233 

0234 

0235 

0236 

0237 

0238 

0235 

C 240 

C 24 1 

0242 

C 243 

0244 

0245 

C 2 46 

C 247 

C 248 

0249 

025C 

0251 

0252 

C 2 5 3 

C 254 

0255 



EDLD * -EOLO 



E0S(1 I ) / 8 . ) 



EOLC • E0LD/10. 

IFIEOSI ll.GT.O.G) 

FLAG - 0*000001 
QZ - A0 S( 1 TYS(KPTS) - TYSI1U/8.) 
QZ - QZ *1000000* 

012 QZ = QZ / 10 * 

FLAG * FLAG*10. 

1 F I QZ - 10.1 813,812,012 
813 IOZ * OZ ♦ 1. 

ZSCAL =■ I QZ 
ZSCAL = ZSCAL*FLAG 
FLAG = 0.000001 
QZ = A0 S I IEDSIKPTS) 

QZ = QZ*1000000. 

014 QZ =« QZ/1G. 

FLAG * F LAG* 10 • 

I F ( QZ - 10. I 022,014,014 

022 102 - QZ ♦ 1. 

YSCAL ■ I QZ 
YSCAL ■ Y$CAL*FLAG 
00 023 1 - l.KPTS 

TYSl I ) - TYSlI) ♦ TLO 

023 EOS « I ) - EDSI1) ♦ EOLO 
wfi - KPTS 

IFIMR.LE.3Q) GO TO 00 
IFIMR.LE. 59 ) GO TO 76 
MQ * MR/20 
MPLOT * MR/M Q 
DO 75 < - 1. MPLOT 

EDSm - EOS { I *M0 1 

75 TYSU) ■ TYS(1*M0) 

MR x MPLOT 
GO TO 80 

76 MD = MR - 30 
MO P 1 * MO ♦ 1 
“ "" 1 



DD 78 
EDS! 1 > 

70 TYSU! 

00 79 
EDS( 1 ) 

79 TYSU! 

MR x 30 
00 CONTINUE 

LA 86 = LA8ELIKC) 



= 1 1 MD 
EOS I l *2 ) 

= TVS ( 1*2) 

I = MOP 1 ,30 
= EDSIUMD) 
TYS< I*MD) 



MAKE A PLOT OF THE OATA USED TO DETERMINE C ANO 0. 

CALL 0RAW<MR,TYS,E0S,1,1,LA8A,KT<TLE,ZSCAL,YSCAL,1, 1,2,2,°, 9,1, LAS 
1 TC I 

PRINT 815, LASTC 

815 FORMAT!//, 1QX, 'LASTC = ’,12) 

TLCW = TLO 
EOLOW = EOLO 
00 824 I = 1 , KPTS 
IUI II = TL NZ ( 1 I ♦ T LOW 
024 YU I I I * FEO< I > ♦ EOLOW 

SUPFR I MPDSE ON THE OATA THE FITTED CURVE FOUND FROM 
SU8RDUT I NE LSGPOL. 

CALL ORAW(KPTS*TU,YU,3»Q»LAS8 ,KT ITLE, ZSCAL, YSCAL , 1, 1, 2,2, 9, 9, 1,LAS 
1 TO I 

PRINT 816, LASTO 

016 FORMAT!//, 10X, 'LASTO = ', 121 
PRINT 817, TLO, EOLO 

817 FORMAT!///, 10X. 'FDR THF PLOT USED TO OETERM! NE"C" AND M 0", THE SCA 
1LES MUST 0 E CORRECTED AS FOLLOWS :',// ,3QX ,' SU8TR ACT ', F4. 2, ' FROM x- 
2 SCALE VALUES', /,30X, 'SUBTRACT', F4. 2, 'FROM Y-SCALE VALUES',//) 

IF THE PREVIOUS GRAPHS HAVE 0E FN FXECUTEO, 00 NOT DRAW 
CREEP CURVES 8UT GO TO STATEMENT 819 ANO PRINT OUT THE 
RESULTS OF THE EMPIRICAL CREEP FIT. 



GO TO 819 
818 CONTINUE 



THE NEXT THREE STATEMENTS PROVIDE THE PROPER VALUE OF THE 
CALLING ARGUMENTS "JC" AND "MODCUR" SO THE CPEEP CURVES 
MAY 0E PLOTTED. 

IFIJC.GT.5I JC = JC - 5 
IFIL.EQ.l) MODCUR = 1 
IFIL.GT.1I MODCUR = 2 
00 66 I = 1 * KR 
66 EMI 1 ) = - El I I 

CALL DR AW < KR, EM, T A, MODCUR, JC.LA8A, IT I TL E , EPS K AL , T SK AL , C , 8 , 2 , 2 , 8 ,15 
lil.LASTA) 

PRINT 29, LASTA 

29 FORMAT!//, 10X, 'LASTA = », 121 

IFIL.EQ.IXI MOOCUR = 3 
IFIL.LT.IX) MODCUR = 2 

LA 88 = LA8ELIKC) 

CALL DRAW I 200, YGM,TG, MOOCUR, 0 ,L ARB, I T I TLE , EP SKAL , TSK AL , 0, 

15,1 , LASTS I 
PRINT 30, LASTB 

30 FORMATI//, 10X, 'LASTB = », 1 2,////) 

819 CONTINUE 



8,2,2, 0, 1 



FITTFD PARAMETERS ARE PRINTED AS OUTPUT. 



THF PARAMFTERS OF THE FITTED CURVE ARE :',//) 
1 2X , ' A3 ' , 12X , 'C * , 13X , * O' ,/J 



PRINT 27 

27 FORMAT I 5X 
PRINT 33 

33 FDRMATI9X, 'Al* ,12X, 'A2 ' ,i^a, ■ 

PRINT 20, A 1 * A2 , A 3 , CC , 0 

28 FORMAT I 1 H ,3014.4, 2F14.4, /// ) 

AA1 = A 1 

A A2 = A 2 
AA3 * A3 
CCC = CC 
QOO = Q 

FITTED PARAMETERS ARE PUNCHED ON CAROS FOR USE IN 
SUBSEQUENT PROGRAMS. 



PUNCH 992, AA1 , AA2, AA3,CCC,QQQ 
992 FORMAT I 3El5.5,2F 10.5) 

00 560 I = 1.NPTS1 

560 TSI I ) = TTI I ) 

THE STATEMENTS ODWN TO STATEMENT 564 ARRANGE FOR A 
SYMBOL "P" TO 9E PRINTED IN THF TABULATED OATA IF AN 
EXPERIMENTAL DATA POINT WAS PLOTTED ON THE CREEP CURVES. 

00 564 I = 1.NPTS1 
J = 1 

561 VAL = A8SITSI I) - TA(JI > 

IFIVAL - 0.000C1) 562,562,563 

562 PMARKH ) = PMKP 
GO TO 564 

563 PMARKH) x ZILCH 
J = J ♦ 1 

IFIJ.GT.KR) GO TO 564 
GO TO 561 

564 CONTINUE 



00 51 I = 1.NPTS1 

YG1II) = E0N(A1.A2,A3 1 CC,0»TT(I ) ) 

51 ERR I I 1 = IIYG1II) - E P I I I )*1QC. ) /EP! I ) 
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0256 

0257 



0256 
0259 
C 260 
0261 
0262 



PRINT OUT THE O0SERVEO STRAIN, FITTED STRAIN, AND 
PERCENT ERROR IN A POINT-8Y -POINT TABULATION. 



u ainu Ktuiitu rim kuin ',A5,» IS LISTED 8 
TIME, HOURS ORSERVFD STRAIN CALC 

"P" IF PLOTTED* ,///) 

= I , NPT S 1 ) 



PRINT 31. RUN 

31 FORMAT! 10X, • THE OATA FITTED ANO PLOTTEO FOR RUN 

IE LOW * ,///»6X, ’POINT NO. 

2 STRAIN rCT ERROR -r- ir ri_ u i i c: u • • t r f 

PRINT 32, < I ,TT(I> , EP( I I.YGKI) , ERR (I) ,PMARK<I) . I 

32 FORMAT! 1H0,I 10 , F 10. 5 , 7X, f 10 . 5 ,<9X , F10. 5 , 5X ,F 8 . 3 , A 12 > 

48 CONTINUE 

42 CONTINUE 
41 CONTINUE 

THE NEXT E OUR STATEMENTS CAUSE ANY "LEFTOVER" DATA $FTS 
0 E FITTED ANO PLOTTEO BEFORE ENDING THE PROGRAM. 

AT THIS POINT THERE ARE BETWEEN ZERO AND (IX - 1) DATA 
SETS NOT YET ANALYZEO. 



0263 


K, 


I F ( NOA T - NOO) 53,53,52 


0264 


52 


IX = NDAT - NDO 


0265 




NDO = NOAT 


0266 




GO TO 40 


0267 


53 


CONTINUE 


0268 




ENO 



fortran 


IV G LEVEL 


1, MOO 1 EON 


OOCi 

0CC2 

0CC3 




FUNCTION EON( Al, A2 , A3,CC. 0. 
IMPLICIT RE AL* 8 (A-Hj. REAL 
IF(Z.EO.0.D> GO TO 704 


0004 




PWR = CC*(Z**0) 


0CC5 




I F ( PWR - 50.) 702,701,701 


CCC6 


701 


TERM = 1.0 


0CC7 




GO TO 7D3 


OCC0 


702 


TERM = 1. - OEXP(-PWR) 


0 CC9 


703 


EON = Al ♦ A 2*Z ♦ A3*TERM 


OOIC 




GO TO 700 


OCli 


7 04 


EON = Al 


OC 12 


7C0 


CONTINUE 


CC13 




RETURN 


OC 14 




END 



08/32/39 



PAGE 0001 



OUTPUT FROM PROGRAM "FITDATA" 



a 2 extrapolation, from subroutine lsqpol 



COEFFICIENTS OF Y=3I + 82*X«-ETC AND ERRORS 
6! 1) = 1. 444382 7D-D2 ERRH= 4.C65E-C4 rt( 2> = 

SIGMA= 1.0256402E-C4 F LEVEL= 0.0 
CCtFFICIENTS OF Y=T 1 *P I *T2*P2+E TC ANO ERRORS 

T! 11= 1.67642860-02 ERRT= 3. 877E-C5 T( 2> = 2. 9625 5530 -05 ER RT= 5.166E-C6 T( 



2. 9625553D-05 ERR8= 5.166E-06 B< 
SUM SO DEV= 7.51384 38E -0 6 



I 



X ( I ) 

7. C0000D0D C 
7. 1 9667C00 0 
7.3833300D 0 
7.78670000 0 
7.80167000 0 
0 . 26 COOOOO 0. 
S.40COCCCO 01 



F (I ) 

1.64000000-02 
1.65500ODD-C2 
1.66D00 000-02 
1 .60500000-02 

1.69000000- 02 

1. 6900000 0- C2 
1.7150000D-02 



Y ( I ) 

1.65176160-02 
1.65750810-02 
1.66311800-02 
l. 6750680D-02 
1.67551150-02 
1.68908980-02 
1.72286290-02 



OEL Y ( I > 

1.17 61616E-04 
2.5880727E-05 
3.1 179792E-05 
-9. 9319 60 IE- 05 
- 1 • 44 8 04 66 E - 04 
-9.1018655E-06 
7 • 86 29440 E-0 5 



W< I ) 

1.0000000E OC 
1 . OOODOGOF 00 
1 • OOOOOCOE OC 
1 . OOOOOOOE DC 
1. OOOOOOOE 00 
1. OOOOOOOE 00 
l.ODDOOOOE OC 
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DETERMINATION OF PARAMETERS 



"C" Ah/D "Q" 



CCtE F 1C I ENTS DE Y*8 1 *82*X*ETC AND ERRORS 
b< 1»* -3.12931050-01 ERRb = 2.330E-02 8( 21= 

SIGMA* I . A 2001 3 3t-01 F LEVEL* 0.0 
CDEEFICIENTS OE Y=T1*PUT2*P2*ETC AND errors 

TI 1>* - 1.A2793AA0-01 E RRT= 2.30AE-02 TI 21* 3.98679070-01 ERRT = 8.2S2F-03 TI 



3.98679070-01 ERR8= 6.252F-03 R< 
SUM SO OE V* 1.9103099E-P2 



1 : 



3 

A 

5 

6 

7 

8 
9 

10 

11 

12 



U : 



X< I ) 

A.09A1AA60 00 
3.68887950 00 
3.A01297AD 00 
3.17821380 OC 
-2.99573230 00 
-2 • 8A 1 81020 00 
-2.70800020 00 
-2.A8A9A67D 00 
-2.30258510 00 
-2.1A8A0580 00 
-2.01A9280D 00 
1.89712000 00 
1 • 69 6A6 750 00 
1.3217A33D 00 
1.0A982210 00 
8. 36255720-01 
6.60350910-01 
2.65698820-01 
1.65325810-02 
A. 16536920-01 
7.01A6251D-01 
1. 10A16350 00 
1.390A6070 00 
1.61277230 OC 
1 • 79 A5389G 00 

2. A9A56A70 00 

3. 15771080 OC 
3.27777A7D 00 
3.39170920 00 
3.82173280 00 
3.87258960 00 
3.91235690 00 
3.951 56A 80 00 
3.9892933D 00 
A.07073A7D 00 
A • 2A8A9520 00 
A. 27620350 00 
A. 30180980 00 



- 1.80 2 3 1 AAD 00 
- 1 • 660A0580 00 
-1 • 5638125D 00 
-1.A7A07730 00 
- 1 . AA559500 00 
-1.39028630 00 
-1.36353890 00 
-1.28620760 00 
-1.18930850 00 
-1.16625190 00 
-1.12098360 00 
-1.055A560D 00 
-9.92617350-01 
-8.55536100-01 
-7.6A026->90-01 
-7. 29183090-C1 
-6.60A62630— 01 
-5. 1A0830 30-01 
-A.06A9696D-C1 
-2.899A25AO-01 
-1.9126AA20-01 
1.62578570-02 
7 . 9A 881 65D-0 2 
1 • 86659A90-0 1 
2.96263960-01 
5. AO 9 36280-01 
6.90561 770-C 1 
8.A2902A7D-C1 
9. 1861O75D-01 
1. 1166AA30 00 
1.1 78660AO 00 
1.26375620 00 
1.37069550 00 
1.A32012AD 00 
1.56382370 00 
1 . A 793 77 50 00 
1.77565370 CO 
1.7A359600 CO 



-1.9A51808D DO 
. 78361010 DO 
.66895710 00 
•580018A0 00 
.50726680 00 
, AA 5901 30 00 
.39255A0D 00 
.30362730 00 
.23092350 00 
• 169A5550 00 
• 1 162A070 00 
.0692731D 00 
9.89277120-01 
8.39882A50-01 
-7.31A7315D-01 
-6. A6328700-01 
-5.76199130-01 
-A. 18859600-01 
-3.06339850-01 
-1 . A6 866500-01 
-3.32726280-02 
1.2727 583D-0 
2.A1A1652D-0 
3. 300A7530-0 
A.0251A060-0 
6.816076AD-0 
9 • A598215D-0 
9. 938A9 120-0 
1.039272A0 00 
1.21071380 00 
1 . 230989A0 00 
1.2A68A380 00 
1.262A751D 00 
1.27751670 00 
1.30998570 CO 
1.38085510 00 
1.39190180 00 
1.A021105D 00 



:t 



OELYI 1) 

1.A286631E-01 
' .2320A23E-01 
_ .051AA62E-01 
1.059A100E-01 
6.1671790E-02 
5.561 A956E-02 
2.9015 l 87E-02 
1 . 7A 19655E-02 
A.161A99AE-02 
3.20 35 A66E -03 
A. 7A289A3 E-0 3 
1.381709AE-02 
3.5A02330E-C3 
1 • 5653651 E-02 
3.25538A0E-02 
8 . 285 A 390 E-02 

8. A263AAAE-02 
9.5223 A27E-02 
1 • 00 1 5708E-0 1 
1 • A307600 E -01 
l • 5799 1 77 E-0 1 
1.1 101 794 E-0 1 
1.6192830F-01 
1.A338803E-01 
1 . 06 25005E-0 1 
1.A067131 E-0 1 
2 • 55A2033E-01 
1.509A662E-01 
I . 2066 1 62 E-01 

9 . AO 6 9481 E— 02 
5 . 232898 1 E-C2 
1.6912 397 E-02 
1.082203AE-01 
1 • 5A A9566E-01 
2.5383794E-01 
9.8522A25E-02 
3.83751 8 7E-0 1 
3 . A 1 A8 550E-0 1 



Will 

l.OOCOOOOE 00 
1 « 000000 OE 00 
1.0000000*= 00 
l.OOCOOOOE OC 
l.OOCOOOOE OC 
l.OOOOOCOE OC 
l.OOOOOOOE 00 
1 . 000000 OE 00 
l.OOCOOOOE 00 
l.OOOOOOOE 00 
l.OOCOOOOE 00 
l.OOOOOOOE 00 
l . OOOOCOOE 00 
1 • OOCOOCOE 00 
l.OOOOOOOE 00 

l.ooenoocE oo 

l.OOOOOOOE 00 
l.OOOOOOOE OC 
l.OOOOOOOE 00 
l.OOCOOOOF 00 
l.OCOOOOOE OC 
l.OOOOOOOE 00 
1. OOOOCOOE OC 
l. OOOOCOOE OC 
l.OOCOOOOE 00 
1 • COOOOC OE 00 
l.OOCOOCCE OC 
l.OOOOOOOE 00 
l.OOOOOOOE 00 
l.OOCOOOOE 00 
l.OOOOOCOE 00 
UCOroPCCE 00 
l.COCOOOOE 00 
l.COOOOCOE OC 
l.OOOOOCOE 00 
1. OOCOOCOE OC 
l.OOOOOCOE OC 
l.OOCOOOOE OC 



THE OATA USEO TO E1N0 "C" AND "0" EOR RUN 2AP13 IS LISTEO 8ELDW: 



. PT. 


LN TIME 


LN I-LNI 1 . — E 0 1 ) 


E ITTEO VALUES 


PCT.ERI 


11 


-A.09A1 


-1.8023 




L.9A52 


-7.9? 


12 


-3.6889 


-1.660A 


— 


1.7836 


-7.A2 


13 


-3.A013 


-1.5638 


- 


1.6690 


-6.72 


1A 


-3.1782 


-1.A7A1 


- 


L.5 800 


-7.19 


15 


-2.9957 


-1.AA56 


— 


L . 5073 


-A. 27 


16 


-2.8A18 


-1.3903 


- 


L.AA59 


-A. 00 


17 


-2.7080 


-1.3635 


— 


1.3926 


-2.13 


18 


-2.A8A9 


-1.2862 


— 


l « 30 36 


-1.35 


19 


-2.3026 


-1.1893 


- 


1.2309 


-3.50 


20 


- 2 • 1 A 8 A 


-1.1663 


- 


L • 16 95 


-0.27 


21 


-2.01A9 


-1.1210 


— 


L . 1 1 62 


0.A2 


22 


-1.8971 


-1.0555 


— ' 


• 0693 


-1.31 


23 


-1.6965 


-0.9928 


-( 


) . 9893 


0.36 


2 A 


-1.3217 


-0.8555 


-( 


3.8399 


1.83 


25 


-1.0A98 


-0. 76A0 


-0.7315 


A. 26 


26 


-0.8363 


-0.7292 


-( 


3. 6A63 


11.36 


27 


-0.660A 


-0.6605 


-( 


3.5762 


12.76 


28 


-0.2657 


-0.51A1 


-( 


3. A 189 


18.52 


29 


0.0165 


-0.A065 


-( 


3.3063 


2A.6A 


30 


0.A165 


-0.2899 


-C 


3. 1A69 


A9.36 


31 


0.7015 


-0.1913 


-( 


3.C333 


82. 6C 



If 


kJ955 


6:8151 




=581: ft 


3 A 


1.6128 


0.1867 


C • 3300 


-76.82 


35 


1.79A5 


C. 2963 


0.A025 


-35.86 


36 


2.A9A6 


0.5A09 


C.6816 


-26.01 


37 


3.1577 


0.6906 


0.9A6O 


-36.99 


38 


3.2778 


0.8A29 


0.9938 


-17.91 


39 


3.3917 


0.9186 


1.0393 


-13. 1A 


AO 


3.8217 


1.1166 


1.2107 


-8.A2 


A 1 


3.8726 


1.1787 


1.2310 


-A.AA 


A 2 


3.912A 


1.2638 


1.2A68 


1.3A 


A3 


3.9516 


1.3707 


1.2625 


7.90 


AA 


3.9893 


1.A32C 


1.2775 


10.79 


A 5 


A.07C7 


1.5638 


1.3100 


16.23 


A6 


A.2A85 


1.A79A 


1.3809 


6.66 


A 7 


A. 2762 


1.7757 


1.3919 


21.61 


A 8 


A. 3018 


1.7A36 


1.A021 


19.59 
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TABULATION OF CALCULATED PARAMETERS AND 
RESULTING CURVE FIJ 



THE PARAMETERS 


OF THE FITTEO CURVE ARE: 








A 1 

0.49570-C2 


A 2 

0. 2963D-04 


A3 C 

0.94870-02 0 


0 

.7313 0. 


3987 




THE OATA 
POINT NO. 


FITTEO ANO 
time, HOURS 


PLOTTED FOR RUN 2AP13 
08SERVE0 STRAIN 


IS LISTEO 8ELOW 
CALC STRAIN 


PCT ERROR 


"P" IF PLOTTED 


! 


o 

o 


0.00496 


0.00496 


C.O 


P 


2 


0.00056 


0.00550 


0.00530 


-3.611 




3 


0.00139 


0.00560 


0.00545 


-2.71R 




A 


0.00278 


0.00575 


0.00560 


-2.633 




5 


0.00417 


0.00590 


0.00571 


-3.277 


P 


6 


0.00556 


0.00600 


0.00579 


-3.440 




7 


0.00694 


0.00610 


0.00587 


-3.B18 




8 


O.OC833 


0.00615 


0.00593 


-3.542 




9 


0.01111 


0.00625 


0.00604 


-3.299 




10 


0.01388 


0.CO63O 


0.00614 


-2.569 




11 


0.01667 


0.00640 


0.00622 


-2.700 




12 


0.025C0 


0.00660 


0.00643 


-2.648 




13 


0.03333 


0.00675 


0.00659 


-2.407 


P 


14 


0. 04166 


0.00690 


0.00672 


-2.546 




15 


C. 05000 


0.00695 


0.00684 


-1.530 





16 


0.05832 


0.00705 


0.00695 


-1.423 




17 


0.06667 


0.00710 


0.00705 


-0.760 




18 


0.08333 


0.00725 


0.00722 


-0.475 




19 


0. 100CO 


0.00745 


0.00736 


-1.173 




20 


0. 11667 


0.00750 


0.00749 


-0.092 




21 


C. 13333 


0.00760 


0.00761 


0.140 


P 


22 


0.15000 


0.00775 


0.00772 


-0.414 




23 


0. 18333 


0.00790 


0.00791 


0.109 




24 


0. 26667 


0.00825 


0.00829 


0.502 




25 


C.35CC0 


0.00850 


0.00859 


1.071 


P 


26 


0.43333 


0. GG 860 


0.00884 


2.779 




27 


0.51667 


0.00880 


0.00905 


2.8*5 




28 


C. 7 66 6 7 


0.C0925 


0.00955 


3.271 




29 


1.01667 


0. CO 960 


0.00993 


3.440 


P 


30 


l. 5167C 


0.01000 


0.01049 


4. 8R0 




31 


2.01670 


0.C1035 


0.01090 


5.290 


P 


32 


3.01670 


0.01110 


0.01149 


3.479 


P 


33 


4.01670 


0.01135 


0.01191 


4.906 




34 


5.01670 


0.01175 


0.C1223 


4. 102 


P 


35 


6.C1670 


0.01215 


0.01250 


2.847 




36 


12.11670 


0.01310 


0.01349 


2. 970 


P 


37 


23.51670 


0.01385 


0.01442 


4.103 


P 


38 


26.51670 


0.01430 


0.01459 


2.048 


P 


39 


25. 7167C 


O.C1455 


0.01476 


1.462 


P 


AO 


45.68330 


C.C1535 


0.01547 


0.758 


P 


41 


48.06670 


0.01550 


0.01556 


0.38C 


P 


42 


50.C167C 


0.01565 


O.G1563 


-0.108 




43 


52.01670 


0.01580 


0.01571 


-0.582 


P 


44 


54 • C 1 6 70 


0.C1590 


C. 01578 


-0.744 


P 


45 


58.60000 


0.C1610 


0.01595 


-0.952 


P 


46 


70.00000 


0.01640 


0.01634 


-0.366 


P 


47 


71.96670 


0.01655 


G. 01641 


-0.870 


P 


48 


73.83330 


0.01660 


0.01647 


-0.794 




49 


77.867C0 


0.01685 


0.01660 


-1.476 


P 


50 


78.01670 


0.01690 


0.01661 


-1 . 739 




51 


82.60000 


0.01690 


0.01676 


-0.855 


P 


52 


84. COG CO 


0.01715 


0.01712 


-0.172 
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X-SCRLE=4 . 00E-03 UNITS INCH. 

V-SCRLE=1.00E+Q1 UNITS INCH. 

HERING, F. S- BOX 21, CREEP OF POLVPROPVLENE 
TIME (HOURS) US. STRAIN ( INCHES/INCH) 
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X-SCPLE=2. 00E-02 UNITS INCH. 
Y-5CPLE=2. 00E+01 UNITS INCH. 



HERING, F- 5- , BOX 21, CREEP OF NYLON 
TIME (HOURS) US. 5TRRIN ( INCHES/INCH) 
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j- 5CHL E=d . 00E-03 UNITS INCH. 
V-SCRL E=1 . 00E+01 UNITS INCH. 



HERING.F-S. BOX 21 , CREEP OF POLYCRRBONRTE 
TIME. HOURS US- STRAIN. INCHES/INCH 
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X-SCRLE=2. OOE+OO UNITS INCH. 
V-SCRLE=6. 00E-01 UNITS INCH. 



HERING.F.S. BOX 21. CREEP OF POLVETHYLENE 
RATE CONST. DETERMINATION LN (-LN (1-ED) ) US- LN TIME 



I 



8 
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APPENDIX C 



COMPUTER PROGRAM FOR PREDICTING CREEP BEHAVIOR OF COMPOSITES 
As discussed in Section 7 s the prediction of composite behavior was 
made by assuming the analytic expression for creep of the components to 
be an equation of state. Input data to the Fortran IV computer program 
called COMCUR consists of the experimental creep data for the composite, 
the volume fractions a and b for each component, and the stress coeffi- 
cients and exponents for each component. These were computed from stress 
correlations of the creep parameters for individual components obtained 
from program FITBATA* For a series of times ranging over the testing 
period, the stress in each component is computed. This is done by a 
Newton-Raphson iteration since the unknown cannot be expressed explic- 
itly as a function of the other variables. With the stress distribution 
determined, the composite creep strain, equal to the strain in each com- 
ponent, is computed using the equation of state assumption. Subroutine 
DRAW is used to display the following output: 

(1) Calculated creep curve of component A at the given stress 

(2) Calculated creep curve of component B at the given stress 

(3) Predicted creep curve of composite AB at the given stress 

for the specified volume fractions of components A and B 

(4) Experimentally observed creep data for composite AB 

The load distribution in the composite during the first 200 seconds of 
the test is also presented in graphical form. Tabulated output and graphs 
for a typical run are given on the following pages. The program COMCUR 
and the related function subprograms are annotated with comment cards. 
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COMPUTER PROGRAM "COMCUR" 



FORTRAN IV G LEVEL 1, MOO 1 

COCJL 
“C2 



GATE = 6BI39 



OCC2 

0003 



0004 
0 005 
0006 
OCC7 
0008 



0GC9 

OO10 



0C1I 

0012 

CC13 

0014 



0015 

0016 
0017 



00 1 B 
CC 19 
0020 
0021 



0022 

0023 

0024 

0025 

0026 



0027 

0028 



0029 

0030 

0031 

0032 
C 03 3 



OC35 
0036 
CC 37 



CC3B 

0039 

CC*C 



0041 

0042 

0043 



CC44 

0C45 

C C 46 

0047 

0048 

0049 

CC5C 

0051 

0052 

0053 

0054 

0055 



HEAL* 8 ITITLEI12I »RUN,KTITLE! 121 
REAL KA .KB.NA.NB.MA .M&. K.M.N, LABEL! 30 I . LABB 

0 1 MENSI ON T 2 ( 200 ) , STRA<206). STRB ( 200 ) , ECOM ( 200 I , E A ( 200 I , 

1) , EME AS ( 2Q0 ) .TMEa£(200)»FRACiA( 2001 »FRACTB!200)jE£alC(200)»ERR(200 

2) ,EMEASM(200) ,E AM! 200) ,EBM! 200) .ECOMM! 200) . TO ( 25 ) . SA ( 25 ) . SB ( 2 5 ) . F A 
3 < i 5 1 ,FB!25 ),sIgG< 200),TG!50I , EG I 50 I ,TY<100),TS( lO0),OA< 100 » ,OBC 100 



ER( 200 
E«R( 200 



R E AO IN TITLES FOR GRAPHICAL OUTPUT 



I * 1,12) 
I » 1,12) 



1,3) 



REAO 5, ( ITITLEI I) , 

READ 5, (KTITLEII), 

5 FORMAT ! 6A8 ) 

REAO 6, (LABEL ( I I , 

6 FORMAT l 3 A 4 ) 

REAO IN THE FOLLOWING OATA: 

RUN = RUN NUMBER OF EXPERIMENTAL COMPOSITE TEST 
NPTS - NUMBER OF CREEP OATA POINTS FOR THE COMPOSITE. 
TMAX = MAX. TIME IN HOURS TO WHICH THE COMPUTED 
CREEP CURVES ARE TO BE ORAWN. 

TSKAL * TIME SCALEIHOURS/ INCH) OF CREEP GRAPHS 
EPSKAL = STRAIN SCALE ( I N/ I N PER INCH) OF CREEP CURVES 

REAO 7, RUN, NPTS, TMAX , TSKAL, EPSKAL, GL 

7 FORMAT ( A 8 • I 12,4F10.0) 

REAO IN EXPERIMENTAL OATA 

EMEAS = DEFORMATION IN MILS 
TMEAS * TIME IN HOURS 

REAO B, (EMEAS(I), I = l, NPTS) 

8 FORMAT! 14F 5. 1 ) 

REAO 9. ( TME A S ( I ) , I = 1 , NPTS) 

9 FORMAT! 7F10. 5) 

REAO IN THE STRESS CORRELATED PARAMETERS OF THE EMPIRICAL 
CREEP EOUATIONS FOR COMPONENTS A AND B 

REAO 10, KA,BA,PA.GA.CA,NA,RA,MA 
REAO 10,KB,BB,PU,0B,CB,NB,RB,MB 

10 FORMAT! 3E15.5.5F6.4I 

REAO IN THE VOLUME FRACTIONS (A ANO B) OF EACH COMPONENT 
ANO THE STRESS LEVEL OF THE EXPERIMENTAL TEST ( STP ) 

REAO 11, A , B 

11 FORMAT! 2F10.0) 

REAO 12, STR 

12 FORMAT(FIO.O) 

GENERATE 200 TIME INTERVALS AT WHICH TO PREDICT THE CREEP 
CURVE. 



TINC = TMAX/ 200 . 
T2< 1) = TINC 

I * 2 , 2C0 

TZ(I-i) ♦ 



00 3- 
31 TZIII 
J = 1 



ESTIMATE THE STRESS IN COMPONENT B AT THE START OF THE 
TEST 

GUESS = STR*B 
00 32 I * 1 » 2CQ 

IN THE FUNCTION SUBPROGRAM "STRESS", COMPUTE THE STRESS 
IN COMPONENT B AT THE 200 TIME INTERVALS !STRB». USE THE 
LAST COMPUTEO VALUE OF STRB AS THE NEW GUESS AT THE NEXT 
TIME INTERVAL. 



»NA,NB, RA*RB,MA»MB, 



STRB! I ) * STRESS(GUE$S,T2 ! I ) , KA , KB , PA , PB , BA , 

1CA , CB , QA , OB , A , R , ST R ) 

IF($TRB4 I I.E0.50000. I GO TO 774 
J * J ♦ 1 

32 GUESS = STRB! I ) 

00 33 I = 1,200 

CALCULATE THE STRESS IN COMPONENT A < STR A ) • 

STRA(I) = ! S TR - B*STRB(I))/A 

USE FUNCTION SUBPROGRAM "STRAIN" TO COMPUTE THE FOLLOWING: 
ECOM = PREOICTEO CREEP STRAIN IN COMPOSITE WHEN 
L CAGE 0 AT STRESS LEVEL EQUAL "STR" 

EA = COMPUTEO STRAIN IN COMPONENT A ALONE WHEN 
LOAOEO AT STRESS EOUAL TO "STR". 

EB = COMPUTEO STRAIN IN COMPONENT B ALONE WHEN LOAOEO 
AT STRESS EOUAL TO "STR". 

ECOM ( 1 ) = STRAIN(STRB( I ),T2( I ),KB,PB,8B,NB,RB,MB.CB,0B) 

EA ( I ) = STKAIN(STR,T2 ( I I , K A, P A , BA , NA , R A , MA , C A , OA ) 

E B ( I ) = STRAIN! STR,T2 ! I ) ,KB,PB,BB,N9,RB,MB,CB,0B ) 

CHANGE THE SIGNS FOR CLARITY ON GRAPHICAL PRESENTATION. 

EAM( I ) = -EA! I ) 

EBM( I ) = -EB! I ) 

33 ECOMM! I ) = -ECCM( I ) 



CHANGE EXPERIMENTAL OEFORMAT ION IN MILS TO STRAIN (IN/IN). 
GL) 



00 34 I = l.NPTS 
EMEAS ! I ) = EMEAS ( I ) / ! 1000 
34 EMEASM! I ) = -EMEAS! I > 



THE NEXT TWELVE STATEMENTS USE SUBROUTINE DRAW TO PLOT THE 
COMPUTEO CREEP CURVE FOR THE INOIVIOUAL COMPONENTS ANO 
FOR THE COMPOSITE. 



0, LABB, I TITLE, EPSKAL, TSKAL, 0,8,2,2,8,15, 1, 

, 12 ,// ) 

2,2,8, 15,1, 



LABB * LABEL! 1) 

CALL ORAW! 200,EAM,T2 
1LASTAI 

PRINT 40 « L ASTA 

40 FORMAT! //, 10X. * LASTA = 

LABb = LABEL ! 2 ) 

CALL OR AW! 200, FBM.TZ, 2. 0,LA8B, I TITLE, EPSKAL, TSKAL, 0 
1LAST8 ) 

PRINT 4 1 , L ASTB 

41 FORMAT! //, 10X. » LASTB * * , I 2 , / / » 

LABB = LABEL!3) 

CALL OR A W( 200, ECOMM, T 2, 2,0, LABB, I T I TLE , EPSKAL , TSKAL ,0 , 8 , 2 , 2 
l.LASTM) 

PRINT 42 , LASTM 

42 FORMAT! // ,10X , » LASTM = * , I 2 , / / 1 



PAGE 0001 
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THE NEXT TWENTY STATEMENTS ARE USEO TO SELECT 30 EXPERI- 

OATA POINTS, THE MAXIMUM AllOWEO IN SU0ROUTI NE 

THE SELECTED EXPERIMENTAL POINTS ARE SUPFRIMPOSEO 



0C56 

0057 

0050 

0059 

0060 
0061 
0062 
0063 
0C66 

0065 

0066 



0076 

D075 



0076 

0077 
0C78 
CO 79 
0000 
0001 
0C82 
0003 

0C66 

D085 

0086 



0103 
0 IC6 
01C5 

01C6 
01C7 
0 1 C 0 
0109 



CllO 
C 1 1 1 



0112 

0113 

0116 

0115 

0116 
C 1 1 7 



C I 18 
0119 



012C 

0121 

0122 
012 3 
0126 
D 1 25 



0126 

0127 



0128 

0129 

C13C 

01J1 



C 132 
0133 



0136 

0135 

0136 

0137 
C 1 30 
0139 
C16C 



MENTAL OATA PO 

ORAW. THE SELI 

ON THE PREOICTED CREEP CURVE. 



IFINPTS.GT.59I 

- ei - N*ir“ 



I PTS 



69 



51 



?°C,1f .*! 
J8 lo 53 



30) 

l , I PTS 
MEASMI I ) 



GO TO 776 
GO TO 50 



EMEA. . r . 

TMEASUl 

. 53 

50 JF • NPTS - 30 
DO 51 I * l.JF 
EGI I I - E ME A SM! 2 *1 ) 

TG!I) - TME AS ( 2*1 ) 

JF2 - ( 2* JF ) ♦ 1 

^,{^1 : ^i N sSIh 

52 TGil-JF) * TMEASlI) 

IPTS • 30 

qnTinuE 

ORAWI I PTS, EG, TG, 3, 3, RUN, IT I TLE , EPSKAL , T SK AL , D , 8 , 2 , 2, 0,15,1 ,L A 

"PRINT 63.LASTD 
63 FORMATI // , lOX, • LASTO * • • I 2 , // I 

GENERATE 1DD TIMES (TY) COVERING THE FIRST 2D0 SECONOS 
OF THE TEST. USE FUNCTION $U0 PROGRAM "STRESS’* TO FIND 
THE CORRESPONDING STRESSES 1 06 1 IN COMPONENT 0. 

FINC = 2./3600. 

TY!1> * FINC 
DO 7D I * 2,100 
70 TY Cl) = TY( 1-1) ♦ FINC 
J * 1 

GUESS = STR*0 
DO 73 I * 1,100 

OB I I 1 = STRE$$!GUESS,TY( I I , KA , K0 , PA , P0 , 0A , 

1C0,OA,O0.A.0,STR) 

IFIC8I I ) .EQ.50D00. ) GO TO 776 
J * J ♦ 1 
73 GUESS = D0U ) 

calculate the stress IN COMPONENT A (OA) ano thf 
LOAD FRACTIONS: 

FRACTA = FRACTION OF THE LOAO CARRIED BY COMPONENT A. 
FRACT0 * FRACTION OF THE LOAD CARRIED 0Y COMPONENT 0. 



,NA,NR,RA,RB,MA,M0 v CA, 



OC07 




DO 76 I = l.IDD 


OC00 




DA! I ) * I STR - B*DB I I ) ) /A 


D089 




FRACTA! I ) * -!DA! I )*A)/STR 


0C90 




FRACTB! I ) * - ! DBI I ) *0 ) /STR 


0C91 




76 TS( I ) * TY! I ) +3600. 


0092 




LABB = LABEL I 1 ) 


0093 




KTITLE! 12) * RUN 




c 


USE SUBROUTINE DRAW TO PLOT 




c 


THE FIRST 200 SECONDS. 


0096 


c 


CALL DRAW! 1DD, FRACTA ,TS,1,O,LA00,KT] 
1STA) 

PRINT 60.LASTA 
LABB = LABELI2) 


0C5 5 




0096 




0C9 7 




CALL DRAW! 100, FRACTB, TS, 3,0, LABB,KT 
1 STB) 


CCS 0 




PRINT 6 1 , L AST0 


0C99 




DO 75 I - 1 , IDO 


0100 




FRACTA! I ) = -FPACTA! I ) 


0101 




75 FRACTB! I ) = -FRACTB! I ) 


0102 




J = 1 



,0,6,2,2,6,15,1 , LA 



,0.2,20. ,0,6,2,2,6,15, 



CALCULATE THE STRESS ISIGG) IN COMPONENT 0 (USING FUNCTION 
SU0PROGRAM "STRESS") CORRESPONDING TO THE EXPERIMENTAL 
TIMES (TMEAS). USE FUNCTION SUBPROGRAM "STRAIN" TO F INO 
THE PREOICTED STRAIN (ECALC) CORRESPONDING TO THE DATA 
POINTS SO A TA0ULATED COMPARISON MAY BE MADE. 

GUESS = STRBIl) 

DO 36 I = 1 , NPT S 

S I GG( I ) * STRESS (GUESS, TME AS I I ) ,KA , K0 , PA , PB, 0A , 0 0 ,NA , N0 , RA , R 0 , M A , 

1 MB ,CA ,CB , QA , QB , A,B * STR ) 

IF(SIGG(I).EQ. 50000.) GO TO 776 

ECALC ( I ) = STRAIN(SIGG( I ),TMEASm,K0,P0,B0,NB,R0,M8,C0,O0) 

J = J + l 

36 GUESS * S IGG ( I ) 

COMPUTE THE PERCENT ERROR IN PREDICTED STRAIN 
DO 37 I = 1 , NPTS 

37 ERR(l) = IIEMEASU) - ECALC! I) ) /E ME AS (I ) ) ♦IDO. 

FINO THE LOAD FRACTIONS IFA,F0) AND THE STRESSES AT 25 
SELECTED TIMES (TO) COVERING THE DURATION OF THE TEST. 



I * 1.25 

= (STRAI I*B-7)*A)/STR 
( S TR0 ( I *0-7 ) *0 ) /STR 



00 30 

FA (I) 

F0( I ) 

TO(I) = TZ I I *8-7 ) 



STRAI I *0-7 ) 

30 S0I I ) = STR0 ( I *B-7 ) 

PRINT OUT THE COMPUTED RESULTS IN TABULAR FORM. 

PRINT 20, RUN 

2D FORMAT! /////, 10X,' COMPARISON OF THE CALCULATED AND EXPERIMENTAL V 
IALUES OF CREEP ST RA I N • , / , 1 OX , » FOR COMPOSITE RUN',A0,' IS GIVEN BE 
2 LOW : ’,///) 

PRINT 21 

21 FORMAT ( 10X , • TIME. HOURS EXPTL. STRAIN CALC. STRAIN 

1 PC T. ERROR*,//) 

PRINT 22, (TMEASI I ) .EMEASI I ),ECALCI I ),ERR(I), I = l.NPTS) 

22 FORMAT! IHO,P2O.5,F10.6,F2O.6,F15.2) 

PRINT 23 

23 FORMAT! /////, 10X,' THE CALCULATED STRESSES CARRIED BY EACH MATERIA 
1L OF THE COMPOSITE AND THE CORR E SPOND I N G * , / , 1 OX , • LOAD FRACTIONS A 
2RE GIVEN BELOW FOR 25 SELECTEO TIMES COVERING THE DURATION OF THE 
3TEST* »/// J 

PRINT 26 

26 FORMAT ! 7X» ' TIME, STRESS IN STRESS IN LOAD FRACT 

1 1 ON LOAD FRACTION' ,/,7X, • HOURS MATERIAL "A" MATERIAL 

2 ”8 " CARRIED BY "A" CARRIED BY "B"',//) 

PRINT 25, (TO! I) , SAID, SB! I ),FA( I ),FB! I ), I = 1,25) 

25 FORMAT! IHO , F 1 3 . 5 ,F 1 5 . I , F 1 7. I , F 16. 3 , F 17 . 3 ) 

PRINT 26 

26 FORMAT! /////, 10X,» THE LOAO FRACTIONS ANO STRESSES CARRIED BY EACH 
1 1 , / , 1DX ■ T COMPONENT DURING THE INITIAL PART OF THE TEST ARE LISTED 
2 BELOW:*,///) 

PRINT 27 

27 FORMAT! 7X. • TIME, STRESS IN STRESS IN LOAD FRACT 

1 1 ON LOAD FRACTION* ,/,7X, ' SECONOS MATERIAL "A" MATERIAL 

2"B" CARRIED BY "A" CARRIED BY "B"'.//) 

PRINT 20, ITS! I ),0A( I ) , UB ( I). FRACTA! I ),FRACT0( I ), I = 1,100) 

28 FORMAT! 1H0,F13.6,F15.1,F17.1,F16.3,F17.3) 

GO TO 776 

776 PRINT 775, J 

775 FORMATI////, 2DX, » CONVERGENCE FAILURE, J =*,I3,///) 

776 CONTINUE 
END 
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PACF 0001 



FORTRAN IV C LEVEL 1, MOO 1 



0CC2 

00C3 

OCCA 

0005 

0CC6 

00C7 

OOCB 

0CC9 

0010 

§811 
0013 
0C14 
OC 15 
OC 16 
OC 17 
CC IB 
0C19 



002C 

00*1 

0C22 



FUNCTION STRESS (SIGB, T,KA,KB,PA,P0, BA , B B , NA , NB-, R A , R R , MA , M0 , C A ,C B , 

THIS SUBPROGRAM USES A NEWTON-RAPHSON ITERATION TO 
COMPUTE THE STRESS IN COMPONENT B I S I GR > FOR A GIVEN 
TIME ( T I • IN THE FUNCTION (FUNCT) THE ONLY UNKNOWN IS 
SIGB. THE FUNCTION IS FORMED BY EQUATING THE STRAIN IN 
EACH COMPONENT. THE STRAIN IN EACH COMPONENT IS GIVEN BY 
THE PREVIOUSLY OETERM 1NEO EMPIRICAL CREEP EOUATION, 

WHOSE PARAMETERS ARE RE AO IN AS ARGUMENTS OF THE 
SUBPROGRAM. THE FUNCTION IS FORMEO BY: 

FUNCT = STRAIN IN A - STRAIN IN B = 0. 

THE DERIVATIVE (OERIV) OF THE FUNCTION IS TAKEN 
ANALYTICALLY. 

1QA.QB,A,0.STR) 

REAL KA»KB,NA,N6»MA,M0 
L = 0 

101 S1GB0L * SIGB 

SIGA = (STR - B*S 1GB ) /A 

TERM1 = KA*ISIGA**NA> - KB* l S 1G0**NB ) ♦ BA* I S I GA **R A » *T - BB*ISIGB 
i**RB» *T 

P WRA = CA*(T**OA» 

IFIPWRA-50.) B2 » 82 » B1 
B1 GO A * 1.0 
GO TO B3 

82 GGA * 1. - EXPl-PWRA) 

B3 PWRB * CB*(T**OB) 

IFIPWRB-50.) B5,B5,B4 
BA GGB = 1.0 
GO TO B6 

85 GGB = 1. - EXP(-PWRB) 

B6 FUNCT = TERM 1 ♦ F A * ( S I GA* * M A ) * GG A - P B* I S I GB* * MB ) *GGB 
H = -B/A 

OERIV = H*(NA*KA*( S1GA**INA-1.I) ♦ R A *BA * ( S I GA** ( RA- 1 . )) * T * ma*PA 
1*(SIGA**(MA-1. )>*GGA) - (NB*KB*ISIG0**(NB-1.) ) ♦ RB*BB* ( S I GB* * ( PB- 
21.»I*T ♦ MB*P0*( SIGB**(M8-1. I)*GGBI 
SIGB = SIGB - I FUNCT /OERIV) 

OIFF = ABSISIGB - SIGBOLI 
L * L *1 

IF THE ROOT OF FUNCT ( = SIGB! IS NOT FOUNO IN TWENTY 
ITERATIONS, ASSUME CONVERGENCE FAILURE. 



CC23 




ifil.gt. 


201 


GO 


TO 


94 


0024 




IFIOIFF 


- 0 


.51 


93, 


i 93,101 


0025 


93 


STRESS = 


S1G0 




0026 




RETURN 










C027 


94 


STRESS = 


50000. 






OC 20 




RETURN 










0 C29 




ENC 











FORTRAN IV G LEVEL 1, MOO 1 



C0C1 

CCC2 



FUNCTION STRAI N( S, 2 
REAL k, p, n 



STRAIN 

K,P,B,N, 



OATE = 60139 



THIS SUBPROGRAM USES THE CREEP EOUATION, WHOSE 
ARt R E AO IN AS SUBPROGRAM ARGUMENTS, TO COMPUTE 
FOR A GIVEN STRESS(S) ANO A GIVEN T f ME I 2 » . 



CCC3 




POWER = C*lZ**0) 


0004 




IF! POWER -50.1 62,62 


OCC 5 


61 


PROO = 1.0 


CCC6 




GO TO 63 


CCC7 


62 


PROO = 1 . — E X P ( -POWER 


OCOB 


63 


STRAIN = K*!S**N> ♦ 1 


CCC9 




RETURN 


0C10 




END 



p*<s**M»*pRon 



PARAMETEPS 
THE STRAIN 



PAGF D001 
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TABULATED OUTPUT FROM 
PROGRAM COMCUR 
FOR PE/PP COMPOSITE 
STRESS = IOOOPSI 



COMPARISON OF 


THE CALCULATEO 


ANO EXPERIMENTAL VALUES OF 


CREEP STRAIN 


FOR COMPOSITE 


RUN AA1 1 IS 


GIVEN BELOW: 


- 


TIME. HOURS 


EXPTL. STRAIN 


1 CALC. STRAIN 


PCT * ERROR 


0.00056 


0.005500 


0.004749 


13.66 



0.00139 


0.006200 


0.005013 


19.15 


0.00270 


0.006400 


0.005274 


17.60 


0.00417 


0.006650 


0.005456 


17.96 


0.00556 


0.006000 


0.005599 


17.66 


0.00694 


0.006900 


0.005719 


17.12 


0.00033 


0.007000 


0.005023 


16.81 


0.01111 


0.007150 


0.006001 


16.08 


0.01380 


0.007300 


0.006148 


15.70 


0.01667 


0.007400 


0.006277 


15.18 


C. 02222 


0.007550 


0.006492 


14.01 


0.02777 


0.007650 


0.006672 


12.79 


0.03330 


0.007750 


0.006826 


11.92 


0.05000 


0.007900 


0.007199 


0.87 


0.00333 


0.000350 


0.007724 


7.50 


0.13333 


0.000650 


0.000264 


4.47 


0.21667 


0.009050 


0.000079 


1.09 


0.30000 


0.009350 


0.009323 


0.29 


0.30333 


0.010050 


0.009674 


3.74 


0.63333 


0.010500 


0.010430 


0.66 


1.00000 


0.011000 


0.0121 22 


-2.73 


3.30000 


0.012900 


0.013128 


-1.77 


3.00000 


0.013050 


0.013361 


-2.30 


4.00000 


0.013500 


0.013741 


-1.79 


6.03330 


0.013900 


0.014108 


-1.49 


10.00300 


0.014600 


0.014904 


-2.09 


21.64999 


0.015900 


0.016035 


-0.05 


23.64999 


0.015900 


0.016165 


-1.67 


27.16699 


0.016250 


0.016371 


i 

o 

-j 


20.79999 


0.016400 


0.016459 


-0.36 


30.04999 


0.016400 


0.016523 


-0.75 


34.29999 


0.016650 


0.016727 


-0.46 


45.63300 


0.017200 


0.017193 


0.04 


40.01669 


0.017300 


0.017281 


0.11 


50.01669 


0.017400 


0.017353 


0.27 


51.06699 


0.017500 


0.017418 


0.47 


53.04999 


0.017550 


0.017487 


0,36 


50.20000 


0.017600 


0.017633 


O' 

0 

1 


69.93300 


0.018000 


0.018003 


-0.02 


71.03299 


0.018000 


0.010061 


-0.34 


73.59999 


0.018050 


0.018114 


-0.36 


75.93300 


0.018100 


0.010184 


-0.46 


77.96700 


0.018200 


0.010243 


-0.24 


61.66699 


0.018250 


0.010351 


-0.55 


93.43300 


0.010550 


0.018685 


-0.73 
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THE CALCULATED STRESSES CARRIED BY EACH MATERIAL OF THE COMPOSITE ANO THF CORR F S POND t MG 
LO AO FRACTIONS ARE GIVEN BELOW FOR 25 SELECTEO TIMES COVERING THE DURATION OF THF TEST 



TIME, 

HOURS 


STRESS IN 
MATERIAL "A" 


STRESS IN 
MATERIAL "8" 


LOAO FRACT ION 
CARRIED BY "A" 


LOAD FRACTION 
CARRIED RY «B» 




C.5000C 


852.6 


1141.0 


0.417 


0.583 




4.50000 


824.9 


1167.5 


0.40 3 


0.597 




8. 50000 


821.2 


1171.1 


0.402 


C.598 




12.50000 


819.9 


1172.4 


0.401 


0.599 




16. 50000 


819.3 


1172.9 


0.401 


C.599 




2C. 50000 


819.1 


1173.1 


0.401 


0.599 




24.50000 


818.9 


1173.3 


0.400 


0.600 




28. 50000 


818.9 


1173.3 


0.400 


0.600 




32.50000 


818.8 


1173.4 


0.400 


0.600 




36.50000 


818.7 


1173.5 v 


0.400 


0.600 


A = PE 


40.50000 


818.6 


1173.6 


0.400 


0.600 


44.50000 


818.5 


1173.7 


0.400 


0.600 


B= PP 


48. 50000 


818.4 


1173.8 


0.400 


C.600 




52.50000 


818.3 


1173.9 


0.400 


0.600 




56.50000 


818.1 


1174.0 


0.400 


0.600 




60.50000 


818.0 


1174.2 


0.400 


0.600 




64.50000 


817.8 


1174.3 


0.400 


0.600 




68.50000 


817.6 


1174.5 


0.400 


0.600 




72.5000C 


817.4 


1174.7 


0.400 


0.600 




76.50000 


817.2 


1174.9 


0.400 


C.6C0 




80.5000C 


817.0 


1175.1 


0.400 


0.600 




84.50000 


816.8 


1175.3 


0.399 


0.601 




88 • 5000 C 


816.6 


1175.5 


0.399 


0.601 




52.50000 


816.4 


1175.7 


0.399 


0.6C1 




56.50000 


816.1 


1176.0 


0.399 


C.6C1 





THE LOAO FRACTIONS ANO STRESSES CARRIED 8Y EACH 
COMPONENT DURING THE INITIAL PART OF THE TEST ARE L 



ISTEO BELOW: 



TIME, 

SECONOS 


STRESS IN 
MATERIAL "A" 


STRESS IN 
MATERIAL "8" 


LOAO FRACTION 
CARRIEO 8 Y "A" 


LOAD FRACTION 
CARRIEO 8 Y "B 


2. 


939.9 


1057.5 


0.460 


0.540 


4. 


936.5 


1060.7 


0.458 


0.542 


6. 


933.9 


1063.2 


0.457 


0.543 


8 . 


931.8 


1065.2 


0.456 


0.544 


10. 


930.0 


1067.0 


0.455 


0.545 


12. 


92 8.4 


1068.5 


0.454 


0.546 


14. 


927.0 


1C69.9 


0.453 


0.547 


16. 


925.7 


1071.1 


0.453 


0.547 


18. 


924.5 


1072.3 


0.452 


0.548 


20. 


923.3 


1073.4 


0.452 


0.548 


22. 


922.3 


1074.3 


0.451 


0.549 


24. 


921.3 


1075.3 


0.451 


C.549 


26. 


920.4 


1076.1 


0.450 


0.550 


28. 


919.6 


1077.0 


0.450 


0.550 


30. 


518.7 


1C77.8 


0.449 


0.551 


32. 


918.0 


1078. 5 


0.449 


0.551 


34. 


917.2 


1079.2 


0.449 


0.551 


36. 


916.5 


1079.9 


0.4^8 


0.552 


38. 


915.8 


1C80.6 


0.448 


C.552 


40. 


915.2 


1C81.2 


0.448 


0.552 


42. 


914.5 


1081. B 


0.44 7 


0.553 


44. 


913.9 


1082.4 


0.447 


0.553 


46. 


913.3 


1083.0 


0.447 


0.553 


48. 


912.7 


1C83.5 


0.446 


0.554 


50. 


912.2 


1084.0 


0.446 


0.554 


52. 


911.7 


1084.5 


0.446 


0.554 


54. 


911.1 


1C85.0 


0.446 


0.5 54 


56. 


910.6 


1085.5 


0.445 


0.555 


58. 


910.1 


1C86.0 


0.445 


C.555 


60. 


909.7 


1086.4 


0.445 


0.555 
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62 . 


909.2 


1086.9 


0.445 


0.555 


64 . 


908.8 


1087.3 


0.444 


0,556 


66 . 


908.3 


1087.7 


0.444 


0.556 


68 . 


907.9 


1088.2 


0.444 


0.556 


TO . 


907.5 


1088.6 


0.444 


0.556 


72 . 


907.0 


1088.9 


0.444 


0.556 


74 . 


904.6 


1089.3 


0.443 


0.557 


76 . 


906.3 


1089.7 


0.443 


0.557 


78 . 


905.9 


1090.1 


0.443 


0.557 


80 . 


905.5 


1090.4 


0.443 


0.557 


82 . 


905.1 


1090.8 


0.443 


0.557 


84 . 


904.8 


1091.1 


0.442 


0.558 


86 . 


904.4 


1091.5 


0.442 


0.558 


88 . 


904.1 


1091.8 


0.442 


0.558 


90 . 


903.7 


1092.1 


0.442 


0.558 


92 . 


903.4 


1092.5 


0.442 


0.558 


94 . 


903.1 


1092.8 


0.442 


0.558 


96 . 


902.7 


1093.1 


0.441 


0.559 


98 . 


902.4 


1093.4 


0.441 


0.559 


100 . 


902.1 


1093.7 


0.441 


0.559 


102 . 


901.8 


1094.0 


0.441 


0.559 


104 . 


901.5 


1094.3 


0.441 


0.559 


106 . 


901.2 


1094.5 


0.441 


0.559 


108 . 


900.9 


1094.8 


0.441 


0.559 


110 . 


900.6 


1095.1 


0.440 


0.560 


112 . 


900.3 


1095.4 


0.440 


0.560 


114 . 


900.1 


1095.6 


0.440 


0.560 


116 . 


899.8 


1095.9 


0.440 


0.560 


118 . 


899.5 


1096.2 


0.440 


0.560 


120 . 


899.2 


1096.4 


0.440 


0.560 


122 . 


899.0 


1096.7 


0.440 


0.560 


124 . 


898.7 


1096.9 


0.439 


0.561 


126 . 


898.5 


1097.2 


0.439 


0.561 


128 . 


898.2 


1097.4 


0.439 


0.561 


130 . 


898.0 


1 C 97.6 


0.439 


0.561 


132 . 


897.7 


1 C 97.9 


0.439 


0.561 


134 . 


897.5 


1 C 98.1 


0.439 


0.561 


136 . 


897.2 


1098.3 


0.439 


0.561 


138 . 


897.0 


1098.6 


0.439 


0 . 561 


140 . 


896.8 


1098.8 


0.439 


0.561 


142 . 


896.5 


1099.0 


0.438 


0.562 


144 . 


896.3 


1099.2 


0.438 


0.562 


146 . 


896.1 


1099.4 


0.438 


0.562 


148 . 


895.9 


1 C 99.7 


0.438 


0.562 


150 . 


895.6 


1099.9 


0.438 


0.562 


152 . 


895.4 


1100.1 


0.438 


0.562 


154 . 


895.2 


1100.3 


0.438 


0.562 



156 . 


855.0 


1100.5 


0.438 


0.562 


158 . 


894.8 


1100.7 


0.438 


0.562 


1 60 • 


894.6 


1100.9 


0.437 


0.563 


162 . 


894.4 


1101.1 


0.437 


0.563 


164 . 


894.2 


1101.3 


0.437 


0.563 


166 . 


894.0 


1101.5 


0.437 


0.563 


168 . 


893.8 


1101.6 


0.437 


0.563 


170 . 


893.6 


1101.8 


0.437 


0.563 


172 . 


893.4 


1102.0 


0.437 


0.56 3 


174 . 


893.2 


1102.2 


0.437 


0.563 


176 . 


893.0 


1102.4 


0.437 


0.563 


178 . 


892.8 


1102.6 


0.437 


0.563 


180 . 


892.6 


1102.7 


0.436 


0.564 


182 . 


892.4 


1102.9 


0.436 


0.564 


184 . 


892.3 


1103.1 


0.436 


C .564 


186 . 


892.1 


1103.3 


0.436 


0.564 


188 . 


891.9 


1103.4 


0.436 


0.564 


190 . 


891.7 


1103.6 


0.436 


0.564 


192 . 


891.5 


1103.8 


0.436 


0.564 


194 . 


891.4 


1104.0 


0.436 


0.564 


196 . 


891.2 


1104.1 


0.436 


0.564 


198 . 


891.0 


1104.3 


0.436 


0 . 564 


200 . 


890.9 


1104.4 


0.436 


0.564 
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X-SCRLE=d . 00E-03 UNITS INCH. 

Y-SCRL E=1 . OOE-tQ J UNITS INCH. 

HERING.B0M 21. COMPARISON OF CRLCULRTED CREEP 
CURUE RND RCTURL DRTR FOR COMPOSITE RB 
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X-SCRLE=2 • 00E-01 UNITS INCH- 
V-SCRLE=2. 00E+01 UNITS INCH. 



HERING, BOX 21, LORD 
OF COMPOSITE RB US- 



FRACTIONS IN EACH COMPONENT 
TIME — CREEP TEST RR11 
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APPENDIX D 



PREDICTION OF STRESS RELAXATION BEHAVIOR FROM EMPIRICAL CREEP EQUATION 

A stress relaxation test is performed by rapidly deforming a material 
to a predetermined strain then observing the decrease in stress as a func- 
tion of time while this strain is maintained. A series of such tests were 
run on the Instron testing machine for the plastics used in this study. 

The empirical equation developed for creep behavior was used to predict 
the relaxation curve. 

(D.i) e= Kcr N + + Pcr^i-x^) 

Again the equation of state assumption was made for predicting relaxation 
behavior. Thus, for a given strain, and a series of times, the corres- 
ponding stress may be found by implicitly solving equation (D.I). This 
was done by a Fortran IV computer program called program RELAX. The input 
data consisted of the creep parameters previously determined and the ini- 
tial strain. A Newton-Raphson iteration was used to solve equation (D.I) 
for stress as a function of time, and these values were tabulated as out- 
put. A plotting subroutine available at the NPGS Computer Facility (SUB- 
ROUTINE DRAW) was used to display the predicted relaxation curves. Fig- 
ure 32 illustrates a computed relaxation curve compared with the experi- 
mental relaxation curve during the initial part of the test. Figure 33 
shows the predicted curve and experimental data over the entire time frame 
of the experiment. These results are typical of all runs; the same general 
trend was observed for each plastic, even at low stress levels. The 
empirical creep equation predicts a greater decrease in stress than that 
observed experimentally. Apparently there is a structure factor influence 
creating irreversibilities that are unaccounted for in the creep equation. 
The equation of state assumption is not valid -- the stress level depends 
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exper\mar»ia\ 
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Predicted Short-Term Relaxation Behavior of Nylon 



Nylon 

= 2.000 psi 




89 



Predicted Long-Term Relaxation Behavior of Nylon Based on Creep Equation 




on the strain history. 

An interesting observation can be made from the experimental data 
in Figure 32. During the first ten seconds of the test, the stress fell 
off rapidly in exponential fashion. This indicates the presence of mole- 
cular mechanisms with relaxation times on the order of several seconds. 

The corresponding retardation time in the creep equation is approximately 
one hour; the predicted curve does not show this sudden decrease at the 
start of the test. It can be concluded, therefore, that the "instantaneous 
strain", parameter in the creep equation, includes a portion of the de- 
layed strain with very short time constants. The time scale of the creep 
test and the method of extrapolation to determine parameter A^ prevents 
the resolution of this rapid but delayed elastic response. From Appendix 
A, Table 3, it is seen that the stress exponent N for the instantaneous 
strain is slightly greater than 1. If the delayed strain component com- 
bined in the term could be removed, this exponent would decrease. 

Although not a necessary consequence of the theory of rubber-like elasti- 
city [2], it appears that the exponent may approach a value of 1, so that 
under instantaneous loading the plastics behave as ideal Hookean solids. 

A modified equation was used with program RELAX, separating the para- 
meter A^ into an instantaneous response and a rapid delayed elastic re- 
sponse : 

w x/ -.n/, „ -t/r*. \ 

(d. 2) € ~ K, 6 + K 2 6 ( l - x ) 

where 

(d.3) K ~ K i + 

The stress coefficients and were estimated from the experiment- 
al relaxation data and the time constant taken to be 10 seconds. 

The predicted relaxation curve then matched the data more closely at short 
times, but the same underestimate of stress level developed at long times. 
(See Figure 33). 
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COMPUTER PROGRAM "RELAX" 



FORTRAN IV G LEVEL 1« MOO 1 



888i 

0003 



0018 

0019 



0C49 

0050 

0051 

0052 

0053 
0C54 
0055 



0056 

0057 
OC 58 
OC 59 
0060 



REAL 
REAL 
DIMFNS 
1200 ) ~ 



K.N.M.LABEL/4H /.LOAD 
*8 I T I TLETi2 ) 

NS I ON PLAST(5).STR<50),T( 
, T8I200) ,TH( 200), AREA I 50) 

i 



OATE * 68125 08/24/38 

/ , LOAD (200) 

200) ,S IG(200)»TA(200)»SIGA(200 ) ,SIG8( 



RE AO IN THE FOLLOWING DATA 

ITITLE ■ TITLE TO BE USEO ON GRAPHICAL OUTPUT 
PLAST » NAME OF MATERIAL TESTEO 

NOAT - NUMBER OF INITIAL STRESS CONOITIONS FOR WHICH 
THE RELAXATION CURVE IS TO BE INVESTIGATED. 

pIrJM * PARAMETERS IN THE EMPIRICAL CREEP EQUATION 
FOR THE MATERIALS TESTEO 
STR = ARRAY OF INITIAL STRESS VALUES 
AREA = ARRAY OE CROSS-SECTIONAL AREAS FOR THE EXPERIMENTAL 
RUNS, CORRESPONDING TO THE INITIAL STRESS VALUES. 



0004 

0005 


13 


RE AO 13.1 ITITLEII ) , I * 1.12) 
FORMAT! 6A8) 


0C06 

0007 


10 


REAO 10 , I PL AST I I ) , I * 1,4), NOAT 
FORMAT I 4A4 , 14 ) 


00C8 

0009 


11 


REAO 11. K.8.P.C,Q.N,R,M 
FORMAT! 3El5.$|5E6. 4) 


0010 

0011 


12 


READ 12 « ( ST R I f ) , I * l.NOAT) 
FORMAT I 8E10. 0) 


0012 

0013 

0014 

0015 

0016 


14 


KCUE = 50 
KCUE * 0 

I F I KCUE • EQ. 50 ) GO TO 16 

lE(KCUE.EQ.O) GO TO 14 

REAO 12 * ( ARE A I I ) , I * l.NOAT) 


0017 


. 16 


CONTI NUE 



THE 00 LOOP "00 102" COMPUTES THE RELAXATION CURVE EROM 
THE EMPIRICAL CREEP EOUATION FOR EACH OE THE INITIAL 
STRESS CONOITIONS. 



00 102 L = 
S * STRIL) 



0020 




EO = 1 


0021 




TINC * 


0022 




Til) > 


0023 




00 32 


0024 


32 


T I I ) : 


0025 




TINC : 


0026 




DO 33 


0027 


33 


T I I ) ; 


0028 




TINC ; 


0029 




00 34 


0030 


34 


Tin 


0031 




TINC : 


0032 




00 35 


0033 


35 


T I I ) : 


0034 




S = 0 



CALCULATE THE INITIAL STRAIN, EO. GENERATE 200 TIMES WITH 
T HE LONGEST ELAPSEO TIME APPROXIMATELY 30 HOURS. STORE 
THESE TIMES IN THE ARRAY "T". 

* I S * *N ) 

1./3600. 

TINC 

I * 2.60 
TII-l) ♦ TINC 
10./3600. 

I * 61,115 
T ( I - 1 ) ♦ T I NC 
1./30. 

I * 116,141 
T I I — 1 ) ♦ TINC 
0.5 

I = 142,199 
TII-1) ♦ TINC 
9 5 *S 

THE DO LOOP "DO 35" USES A NEWTON-R APHSON ITERATION TO 
SOLVE THE EMPIRICAL CREEP EOUATION IMPLICITLY FOR STRESS 
AT EACH TIME IN THE ARRAY "T". THE INITIAL STRAIN IS EO. 
THE FIRST ESTIMATE OF THE RELAXED STRESS AT ELAPSED TIME 
Til) = 1 SECONO IS 0 . 95* S . THE FIRST ESTIMATES AT LATER 
TIMES IS THE RELAXEO STRESS AT THE PREVIOUS ELAPSED TIME. 
STRESSES COMPUTEO FOR THE TIMES T I I ) ARE STORED IN THE 
ARRAY NAMED "SIG". 



0035 




DO 36 I = 1,199 


0036 




LC = 0 


0037 




DIFE = 40000. 


0038 




PWR * C * ( T I I )**0) 


0039 




IEIPWR - 50.) 82,82 


0040 


81 


0 = 1.0 


0041 




GO TO 83 


0042 


82 


0=1.- EXPI-PWR) 


0043 


83 


LC = LC ♦ 1 


0044 




SOLD = S 


0045 




EUNCT = -EO ♦ K*!S*’ 


0046 




OERIV = N*K*(S**IN-: 


0047 




S = S - IEUNCT/DER I' 


0048 




DIEF = ABSISOLD -S ) 



♦ B* (S**R >*T I I ) ♦ P*IS**M)*D 



IE THE STRESS IS NOT EOUNO WITHIN 0.1 PS I IN TWENTY 
ITERATIONS, ASSUME CONVERGENCE FAILURE. 

IEILC.GT.20) GO TO 87 
IE(S.LE.O.O) S * S/2. 

IEIDIEF - 0.1) 36,83,83 

87 PRINT 88. I-LC 

88 EORMATUOx,* CONVERGENCE FAILURE, I = »,I2,‘ LC = NI2,//) 

S = 0. 75*ST R I L ) 

16 SIG(I) = S 

COMPUTE THE "FRACTION OF RELAXATION", EQUAL TO RELAXED 
STRESS OIVIOEO 8Y INITIAL STRESS, AND STORE IN THE 
ARRAY NAMEO "SIGA". 

DO 38 I = 2.199 
TA( ( ) = T(I-l) 

38 SIGAI I ) = SIGI I— 1I/STRIL) 

T A I 1 ) = 0.0 
SIGAI 1) = 1.0 

USE SUBROUTINE ORAW TO PLOT THE PREDICTED RELAXATION CURVE 
OE "SIGA VS. ELAPSEO TIME". 

,1. 1,2, 2, 9, 9,1 .LAST) 



0061 




CALL DRAW! 199. TA,S IGA, 0,0, LABEL, ITITLE, 5 


0062 




PRINT 40, LAST 


OC 6 3 


40 


FORMAT! //,10X.» LAST = ', 12,//) 


0C64 




DO 39 I = 2.200 


0065 




T 8 I I ) = TII-l)*3600. 


0066 




THI I) = T I I — 1 ) 


0067 


19 


SIG8I I ) = SIGI 1-1) 


0068 




SIG8I l ) = STRIL) 


0069 




T8I 1) = 0.0 


007C 




TH(l) = 0.0 


0071 




PRINT 44 


0072 


44 


FORMAT! 1H1) 


0073 




IEIKCUE.E0.50) GO TO 48 



PAGE con 
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0074 

0075 

0076 

0077 



0070 

0C79 

0080 

0C81 

0082 

0003 

0084 



0005 

0086 

0007 

0088 



COMPUTE THE PREQICTEO AELAXEO LOAO BAS EO ON SPECIMEN AREA 
SO THAT EXPERIMENTAL LOAO-TIME OATA MAY BP COMPARE 0 WITH 
PRFDICTEO LOAD-TIME BEHAVIOR. 



00 31 I = 1.200 
31 LOAOI I > = SIGB! I )*AREA( L ) 

PRINT 47.AREAIL) 

47 FORMAT!///, 20X, 'THE SPECIMEN AREA IS* 



F 7 . 5 , * SOUARF INCHES’ ,////) 



PRINT OUT THE PREOICTEO STRESS RELAXATION CURVF IN 
TABULAR FORM. 

PRINT 45, ( PLAST (II. I - 1.4) 

PRINT 49, (TH( I ) , T 8 ( I) , S I GB f I ) , LOAD ( I), I = 1,200) 

49 FORMAT! 1H0.20X, F10. 6,1 3X,F8. 0,1 3X,F 0.1, 13 X,F 0.1 ) 

GO TO 103 
48 CONTINUE 

PRINT 45. I PL AST ! I ) . I - 1,4) 

45 FORMAT! 10X, 'THE CALCULATEO STRESS RELAXATION BEHAVIOR FOR *,4A4.* 
IIS LISTEO BELOW:*. //,20X, 'TIME, HOURS', 10X, 'TIME, SECONDS*, 10X, • STR 
2E SS, PSI ' , 10X . ' LOAD*.//) 

PRINT 46, ( THU ) ,TB( I I ,S IGBI I ) , I = 1.200) 

46 FORMAT! 1HO,2OX,F1O.6,13X,F0.O,13X,F0. 1) 

103 CONTINUE 

102 CONTINUE 
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TABULATED OUTPUT 



FROM "RELAX" 



THE SPECIMEN 


AREA ISO. 04992 SOUARE INCHES 








CALCULATED STRESS 


RELAXATION BEHAVIOR FOR 


NYLON IS LI 


ISTFO 8EL0W: 




TIME, HOURS 


TIME , SECONOS 


STRESS, PSI 


LCAO 




0.0 


0.0 


2000.0 




99. 8 


0.000270 


l. 


1999,1 




99.8 


0.000556 


2. 


1998.3 




99.7 


0.000833 


3. 


1997.5 




99.7 


O.OOliil 


4. 


1996.7 




99.7 


0.001389 


5. 


1996.0 




99.6 


0.001667 


x 6 • 


1995.2 




99.6 


0.001944 


7. 


1994.5 




99.6 


0.002222 


8. 


1993.7 




99.5 


0.002500 


9. 


1992.9 




99.5 


0.002778 


10. 


1992.2 




99.4 


0.003056 


11. 


1991.5 




99.4 


0.003333 


12. 


1990.7 




99.4 


0.003611 


13. 


1990,0 




99.3 


0.003889 


14. 


1989.3 




99.3 


0.004167 


15. 


1988.5 




99. 3 


0.004444 


16. 


1987.8 




99.2 


0.004722 


17. 


1987.1 




99.2 


0.005000 


18. 


1986.4 




99.2 


0.005278 


19. 


1985.7 




99. 1 


0.005556 


20. 


1984.9 




99.1 


0.005033 


21. 


1984.2 




99.0 


0.006111 


22. 


1983.5 




99.0 


0,006389 


23. 


1982.8 




99.0 


0.006667 


24. 


1982.1 




98.9 


0.006944 


25. 


1981.4 




98.9 


0.007222 


26. 


1980.7 




98,9 


0.007500 


27. 


1980.0 




98.8 


0.007778 


28. 


1979.3 




98.8 


0.008056 


29. 


1978.6 




90.8 


0.008333 


30. 


1977.9 




98.7 


0.008611 


31. 


1977.2 




98.7 


0.008889 


32. 


1976.5 




98.7 


0.009167 


33. 


1975.8 




98.6 


0.009444 


34. 


1975.2 




98.6 


0.009722 


35. 


1974.5 




98.6 


0.010000 


36. 


1973.8 




98.5 


0.010278 


37. 


1973.1 




98.5 


0.010556 


38. 


1972.4 




98.5 


0.010833 


39. 


1971.7 




98.4 


0.011111 


40. 


1971.1 




98.4 


0.011389 


41. 


1970.4 




98.4 


0.011667 


42. 


1969.7 




98.3 


0.011944 


43. 


1969.0 




98.3 


0.012222 


44. 


1968.4 




98.3 


0.012500 


45. 


1967.7 




98.2 


0.012778 


46. 


1967.0 




98.2 


0.013056 


47. 


1966.4 




98.2 


0.013333 


48. 


1965.7 




98.1 


0.013611 


49. 


1965.0 




98.1 


0.013889 


50. 


1964.4 




98.1 


0.014167 


51. 


1963.7 




98.0 


0.014444 


52. 


1963.1 




98.0 


0.014722 


53. 


1962.4 




98.0 


0.015000 


54. 


1961.7 




97.9 


0.015278 


55. 


1961.1 




97.9 


0.015555 


56. 


1960.4 




97.9 


0.015833 


57. 


1959. R 




97,8 


0.016111 


58. 


1959.1 




97.0 


0.016389 


59. 


1958.5 




97.0 


0.016667 


60. 


1957.8 




97,7 


0.019444 


70. 


1951.4 




97.4 


0.022222 


80. 


1945.0 




97. 1 



I 
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0.025000 


90. 


1938.8 


96. 8 


0.027778 


100. 


1932.7 


96.5 


0.030555 


110. 


1926.7 


96.2 


0.033333 


120. 


19?0.7 


95.9 


0.036111 


130. 


1914.8 


95.6 


0. 038889 


140. 


1909.1 


95.3 


0 .0*1667 


150. 


1903.3 


95.0 


0.044444 


160. 


1897.7 


94. 7 


0.047222 


170. 


1892.2 


94.4 


0.050000 


180. 


1 886.7 


94. 7 


0.052778 


190. 


1881.2 


93.9 


0.055555 


200. 


1875.9 


93.6 


0.058333 


210. 


1870.6 


93.4 


0.061111 


22 0. 


1865.3 


93.1 


0.063889 


2 30. 


1860.2 


92.9 


0.066667 


240. 


1855.1 


92.6 


0.069444 


'250. 


1850.0 


92.3 


0.072222 


260. 


1845.0 


92.1 


0.075000 


270. 


1840.1 


91.8 


0.077778 


280. 


1835.2 


91.6 


0.080555 


290. 


1830.3 


91.4 


0.083333 


300. 


1825.5 


91.1 


0.086111 


310. 


1820.8 


90.9 


0.088889 


3 20. 


1816.1 


90.7 


0.091666 


330. 


1811.5 


90.4 


0.094444 


340. 


1 806.9 


90.2 


0.097222 


350. 


1802.3 


90.0 


0.100000 


360. 


1797.8 


87, 7 


0.102777 


370. 


1793.4 


89.5 


0.105555 


380. 


1789.0 


89.3 


0.108333 


390. 


1784.6 


89.1 


0.111111 


400. 


1730.2 


38.9 


0.113888 


410. 


1776.0 


88. 6 


0.116666 


420. 


1771.7 


88.4 


0.119444 


430. 


1767.5 


88.2 


0.122222 


440. 


1763. ^ 


8 8.0 


0.124999 


4 50. 


1759.7 


87,8 


0.127777 


460 . 


1755.1 


87.6 


0.130555 


470. 


1751.0 


87.4 


0.133333 


480. 


1747.0 


87.2 


0.136110 


4°0. 


1743.0 


87. n 


0.138888 


500. 


1739.1 


P 6.8 


0.141666 


510. 


1735. 1 


86.6 


0. 144444 


5 20. 


1731.7 


86.4 


0.147221 


530. 


1727.4 


86. 2 


0 . 149099 


540. 


1723.6 


86. 0 


0.152777 


550. 


1719.8 


85.8 


0.155555 


560. 


1716.1 


8^.7 


0.158332 


570. 


1712.3 


85.5 


0.161110 


580. 


1 708.6 


85.3 


0.163888 


590. 


1705.0 


85.1 


0. 166666 


600. 


1701.7 


84. 9 


0.169443 


610. 


1697.7 


84.7 


0.202777 


730. 


1656.7 


82.7 


0.236110 


850. 


1619. 1 


80.8 


C . 269443 


970. 


1584.6 


79. 1 


0.302777 


1090. 


1552.8 


77.5 


0.336110 


1210. 


1573.2 


76.0 


0.369443 


1330. 


1496.7 


74.7 


0.402777 


1450. 


1469.9 


73.4 


0.436110 


1570. 


1445.8 


. 7? « 2 


0.469443 


1690. 


1427.1 


71.0 


0.502777 


1810. 


1401.7 


70.0 


0.536110 


19 30. 


1381.5 


69.0 


0.569443 


2050. 


1362.4 


68.0 


0.602776 


2170. 


1 344.7 


67.1 


0.636110 


2290. 


1327.1 


66.2 


0.669443 


2410. 


mo.-' 


65.4 


0.702776 


2530. 


1295.0 


64.6 


0.736110 


2650. 


1280.1 


63.9 


0.765443 


2770. 


1265.9 


63.2 


0.802776 


2890. 


1252.2 


62.5 


0.836110 


3010. 


1239.2 


61.9 
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0.869443 


3130. 


1226.6 


61.2 


0,902776 


3250. 


1214.6 


> 

o 


0. 936109 


3370. 


1203.0 


60.0 


0.969443 


3490. 


1191.9 


59.5 


1 .002773 


3610. 


1181,2 


*9. 0 


1.036108 


3730. 


1170.8 


58.4 


1.536108 


5530. 


1050.9 


52.5 


2.036108 


7330. 


972.3 


48.5 


2.536108 


9130. 


916.7 


V* 

00 


3.036108 


10930. 


875.3 


43.7 


3.536108 


12730. 


843.5 


42.1 


4.036108 


14530. 


818.4 


40.9 


4.536108 


16330. 


798.2 


39.8 


5.036108 


18130. 


781.7 


39.0 


5.536108 


19930. 


768.2 


38.3 


6.036108 


21730. 


756.8 


37.8 


4.536100 


23530. 


747.3 


37.3 


7.036100 


25330. 


739.3 


36.9 


7.536108 


27130. 


732.4 


?6. 6 


8.036108 


28930. 


726.6 


36.3 


8.536108 


30730. 


721.5 


36.0 


9.036108 


32530. 


717.2 


35.8 


9.536108 


34330. 


713.5 


35.6 


10.036108 


36130. 


710.2 


35.5 


10.536108 


37930. 


707.4 


35.3 


11.036108 


39730. 


704.9 


35.2 


11.536108 


41530. 


702.7 


35.1 


12.036108 


43330. 


700.9 


35.0 


12.536108 


45130. 


699.2 


34.9 


13.036108 


46930. 


697.7 


34. 8 


13.536108 


48730. 


696.4 


34.8 


14.036108 


50530. 


695.2 


34.7 


14.536108 


52330. 


694.2 


34.7 


15.036108 


54130. 


693 . 3 


34.6 


15.536108 


55930. 


692.5 


34.6 


16.036102 


57730. 


691.7 


34.5 


16.536102 


59530. 


691.1 


34.5 


17.036102 


61330. 


690.5 


34.5 


17.536102 


63130. 


689.9 


34.4 


18.036102 


64930. 


689.4 


34.4 


18.536102 


66730. 


689.0 


34.4 


19.036102 


68530. 


688.6 


34.4 


19.536102 


70330. 


688.2 


34.4 


20.036102 


72130. 


68 7,8 


34.3 


20.536102 


73930. 


687.5 


34.3 


21.036102 


75730. 


687.2 


34.3 


21.536102 


77530. 


686.9 


34.3 


22.036102 


79330. 


686.6 


34.3 


22.536102 


8U30. 


686.4 


34.3 


23.036102 


82930. 


686.1 


34.2 


23.536102 


84730. 


685.9 


34.2 


24.036102 


86530. 


685.7 


34.2 


24.536102 


88330. 


685.5 


34.2 


25.036102 


90130. 


685.2 


34.2 


25.536102 


91930. 


685.0 


34.2 


26.036102 


93730. 


684.8 


34.2 


26.536102 


95530. 


684.7 


34.2 


27.036102 


97330. 


684.5 


34.2 


27.536102 


99130. 


684.3 


34.2 


28.036102 


100930. 


684.1 


34.1 


28.536102 


102730. 


683.9 


34.1 


29.036102 


104530. 


683.8 


34.1 


29.536102 


106330. 


683.6 


34.1 


30.036102 


108130. 


683.4 


34.1 



GRAPH TITLED 

HERING.F.S., SOX 21, STRESS RELAXATION OF NVLON 
FRACTION OF INITIAL STRESS VS. TIME, HOURS 
HAS BEEN PLOTTEO. 



LAST * 0 
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APPENDIX E 



PREDICTION OF STRESS-STRAIN CURVES USING THE EMPIRICAL CREEP EQUATION 
A simple method for predicting creep curves from stress-strain curves 
is outlined in Figure 34. As indicated in Figure 34(a), tensile tests are 
run at a series of constant strain rates. A cross-plot. Figure 34(b), is 
constructed for reciprocal strain rate vs. strain at a constant stress. 

By the graphical integration shown in Figure 34(b), times corresponding to 
certain strains may be found for a given stress level, and creep curves 
plotted as illustrated in Figure 34(c). This method implies the existence 
of an equation of state relationship between strain, stress, and strain 
rate, but does not explicitly state such an equation. A similar scheme 
was proposed by Marin and Pao [30] , who obtained good agreement from test 
data on polymethylmethacrylate. 

In this study, the empirical creep equation was used to predict stress- 
strain curves by reversing the procedure outlined in Figure 34. From the 
empirical creep equation, 



The slope (strain rate) at any time for a given stress may be found by: 



For a specified strain rate and a given stress, equation (E.2) may be 



to obtain the corresponding strain. A Fortran IV computer program, called 
SIGEPS, was used to perform these calculations. Output from SIGEPS con- 
sisted of stress-strain data for a material for a given strain rate. 





solved implicitly for time, and this result substituted in equation (E.l) 
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Figure 34 

Prediction of Creep Curves from Stress-Strain Data 
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Program SIGEPS and typical output is given on the following pages* Figure 
35 displays this output along with experimental tensile test data obtained 
at the same strain rate. The same characteristic was observed for most 
materials tested; the predicted tensile strain is greater than the observed 
strain for a given stress level* The deviation from the predicted curve 
becomes greater at higher stress levels* In the tensile test, stress 
changes with time, and internal structure changes probably develop* These 
changes are not accounted for in the creep equation used for the predic- 
tion. In addition, it should be noted that the tensile tests took place 
over a time span of about 10 minutes whereas the time duration of the creep 
test was 100 hours. Consequently the prediction of stress-strain curves is 
based on the creep curve at very short times. As indicated in the tabulat- 
ed output from program FITDATA (See Appendix B) , it is at these short times 
where poorest fit of equation (E.1) usually occurred. 

Figure 36 illustrates an exception to the previous discussion. The 
predicted stress-strain curves for polycarbonate showed excellent agree- 
ment with tensile test data. It is possible that structure factor effects 
in polycarbonate are negligible at the stress levels encountered. Further- 
more, output from FITDATA for polycarbonate showed good agreement between 
the empirical creep equation and the test data at times under 10 minutes. 

It would be useful to further examine the feasibility of this method 
by utilizing an empirical creep equation specifically designed to give 
good data fits at short times. 
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Stress , psi x 10” 




Figure 35 

Predicted Stress-Strain Curve for Nylon 
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Figure 36 

Predicted Stress-Strain Curve for Polycarbonate 
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COMPUTER PROGRAM "SIGEPS 



FORTRAN IV G LEVEL 1, MOO l 



OATE * 68132 



8881 



0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 



0013 

0014 

0015 

0016 



0017 

0018 

0019 

0020 
0021 



0022 

0023 

0024 

0025 

0026 



0027 

0028 

0029 

0030 

0031 

0032 

0033 
00 34 

0035 

0036 

0037 

0038 
00 39 

0040 

0041 

mi 

0044 

0045 

0046 



p£MEN§ION EOO T ( 50 1 , E ( 50 1 * T ( 50) » S IG ( 50 ) , E P ( 50 , 50 1 , PL AST I 10 1 

REAu O. IKifl 
8 FORMAT I4A4 ) 

REAO IN THE FOLLOWING DATA! 

JPTS * NO. OF DIFFERENT STRAIN RATES FOR WHICH STRESS- 
STRAIN CURVES ARF TO 8E COMPUTED 
NOAT = NUMBER OF OIFFERENT STRESSES AT WHICH THE STRAIN 
IS TO BE EVALUATED. 

A1,A2,A3,C,Q - PARAMETERS IN THE EMPIRICAL CREEP 
EQUATION FOR A GIVEN STRESS 
EOOT * ARRAY OF STRAIN RATES FOR WHICH THE STRESS - 
STRAIN CURVE IS TO 8 E DRAWN 
STR =* STRESS CORRESPONDING TO THE CREEP EQUATION 

parameters 

READ 10. JPTS, NOAT 
10 F0RMAT(2I5) 



REAO 14, IE 00 T( II, I = 1,JPTS> 
9 REAO 12, A1,A2,A3.C,Q 
12 FORMAT ( 3E 1 5 • 5, 2F10.0) 



REAO 13, STR 

13 FORMAT ( F10.0 ) 

14 FORMAT! 7F10. 01 

ESTIMATE THE TIME AT WHICH A STRAIN RATE EOOTin OCCURS ON 
THE CREEP CURVE. 

GUESS = !2.*A1)/(E OOT (11*60.) 

00 31 I = 1 , J PT S 

85 PRINT 86, I 

86 FORMAT! // ,20X, * I = ♦, 121 

USE FUNCTION SU8PR0GRA M "TIME" TO FINO THE TIME AT WHICH 
A STRAIN RATE EOOT ( I ) OCCURS FOR THE CREEP CURVE GIVEN BY 
THE EQUATION WHEN STRESS - STR. THEN USE FUNCTION SUB- 
PROGRAM "STRAIN" TO COMPUTE THE STRAIN FOR THE GIVEN 
STRESS. 



T ( 1 1 * TIME! GUESS|A1|A2,A3,C|Q, EOO Till) 
E ( I ) = STRAIN!A1,A2,A3,C,Q,T(I h 
IF(T( I ).EQ. 10000. ) GO TO 771 

- E(I) 

STR 



EP (LC.I ) 
SIG(LC) 



CONTINUE THE PROCESS FOR "NOAT" OIFFERENT STRESS LEVELS. 

31 GUESS = Till 
40 LC = LC+ 1 

I F ( LC - NOAT) 9,9,23 

23 PRINT 24 

24 FORMAT! 1H1 ) 

PRINT OUT THE PREOICTEO STRESS-STRAIN CURVE IN TABULAR 
FORM FOR EACH STRAIN RATE. 

PP INT 25, ! PL A ST ( I ) , I = 1,4) 

25 FORMAT ( 10X, * THE STRESS-STRAIN OATA FOR* ,4A4, * FOR VARIOUS STRAIN 
1ATES IS LISTEO BELOW:',///) 

PRINT 26 

- 26x, , STRAIN rates, in/i n/min» , // ) 



26 FORMAT ( 3X, * STRESS' 

PRINT 27, ( EOOT ( I ) , \ = 1,JPTS) 

27 FORMAT! l2x,5F15. 6,//) 

00 32 J = 1 , NOA T 

PRINT 28, SIG! J),! EP(J, I) , I * 
32 PRINT 29 



1 , JPTS ) 



28 FORMAT! ///.F10. 0,5F15. 6) 

29 FORMAT! 1H0) 

GO TO 42 

771 II = I 

PRINT 772, LC , I 

772 FORMAT! ///, 10X, » CONVERGENCE FAILURE AT LC 

773 gg,L 7 £?M , .V.S"’ PIS 
GO TO 40 

42 CONTINUE 
ENO 



PAGF onoi 
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FORTRAN IV 
OOOI 



LEVEL I« MOO 1 STRAIN OATE = 68132 Ofl/I 
FUNCTION STRAINIA1 ,A?,A3,C,0,Z) 

C FOR A GIVEN TIME <Z) THIS FUNCTION SUBPROGRAM USES THE 
C EMPIRICAL CREEP EOUATION FOR A GIVEN STRESS TO COMPUTE 
C THE CORRESPONDING VALUE OF STRAIN. 



0002 




POWER = 


C*IZ**0) 


0003 




IF! POWER 


- 50.) 


0004 


61 


PROOT = 


1.0 


0005 




GO TO 63 


I .-EXPI-l 


0006 


62 


PROOT = 


0007 


63 


STRAIN = 


A 1 ♦ A2 


0008 




RETURN 




0009 




ENO 





62 » 62 » 61 



P AO F 000 



FORTRAN IV G LEVEL I, MOO 1 



OATE = 6B I 3? 



0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 



00 1 R 

0019 

0020 
0021 
0022 

0023 

0024 



FUNCTION TIME! X, A 1 , A2 , A3 , C , 0 » E 01 



THIS FUNCTION SUBPROGRAM COMPUTES THE T l MF IX) AT WHICH A 
GIVEN STRAIN RATE IEO) OCCURS, 8ASE0 ON THE CREEP PARAM- 
ETERS Al f A2.A3,C.O FOR A GIVEN STRESS. THE FUNCTION 
CALLEO "EUNCT" IS THE DERIVATIVE OF THE CREEP EOUATION, 
OE/OT. MINUS THE GIVEN STRAIN RATE»EO. THF SUBPROGRAM 
USES A NEWTON-RAPHSON ITERATION TO E I NO THE UNKNOWN TIME 
(X), WHICH IS THE ROOT OF THE EOUAT I ON M FIJNCT = 0". 

"OERIV" IS USE 0 IN THE N/R ITERATION AND IS THE OERIVATIVF 
OF "FUNCT" TAKEN ANALYTICALLY. 

L = 0 

R = Q - 1. 

90 XOLO = X 

PRINT R7.L ,XOLO 

87 E0RMATI25X.* L = , ,I2,» XOLO = *,E15.9) 

PWR = C*I X**0) 

IF ( PWR - 50.) 80,80,82 

80 PROO = A3*EXP< -PWR ) 

EUNCT = A2 - EO + ( C *0*PR00* < X ♦* R ) > 



OERIV = PR00*C*0*IR*< X**IR-1. ) ) 
X = X - IFUNCT/OER IV) 
IEIX.LT.O.O) X =X0LD/2. 

OIFF = A8SIX - XOLO) 

PRINT 93, OIFF, FUNCT, OERIV, PWR 
93 FORMAT ( 25 X , * OIFF = »,F15.9,» 
1,» PWR =• '.F10.6J 
L = L ♦ 1 



- C*0*(X**(2.*R) ) ) 



, E 1 5 . 9 , 1 OERIV = 



IF more THAN 20 ITERATIONS ARE USEO TO FINO THE 
TIME IX), ASSUME CONVERGENCE FAILURF. 



IEIL.GT.20) GO TO 82 
IFIOIFF - .00001) 81,81,90 
R1 TIME = X 
return 

82 TIME = 10000. 

RETURN 

ENO 



TABULATED OUTPUT FROM "SIGEPS" 



THE STRESS-STRAIN OATA FOR POLYCARBONATE FOR VARIOUS STRAIN RATES IS LISTEO BELOW: 
STRESS STRAIN RATES, IN/IN/MIN 



0.013040 

2000. 0 . 006QB 7 



2500* 0.007790 



3000. 0.009535 



3500. 0.011291 



4000. 0.013212 



4 500 . 0.015322 



5000. 0.01H104 
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APPENDIX F 



THE BOX DISTRIBUTION FUNCTION FOR RETARDATION SPECTRA 
The use of a continuous distribution of retardation times to character- 
ize delayed elastic response was discussed in Section 6. Not only do 
theoretical considerations favor this approach, but mathematical solutions 
to viscoelastic problems may be considerably simplified, particularly if 
the functions representing the spectra can be handled by transform techni- 
ques. As an example, an analytic expression for creep behavior will be 
derived here, using a box distribution of retardation times as illustrated 
in Figure 37. 

An expression for creep strain corresponding to the four-element model 
of Figure 10, page 25 , may be written: 

(f . i) e = A i + Ait + A 3 (i - JL bt ) 

where the constant b is the rate constant equal to the reciprocal of the 
characteristic relaxation time / ~T , and coefficients A^, A^, and A^, de- 
pend on the material and the applied stress. For a material with many dis- 
creet retardation times, the expression for creep strain given by equation 
(F.l) becomes: 

(f • 2 ) g = A, + Ait + A 

where the constant bl is the rate constant associated with the 

retardation time '71 , and the factors CL i are weighting constants in- 

* 

dicating the relative contribution of the L mechanism to the delayed 
strain. Since the factors CL i are chosen so that — 1 > equation 

(F.2) may be rearranged as: 
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Figure 37 

Box Distribution of Rate Constants 




Figure 38 

Wedge Distribution of Rate Constants 
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(f . 3) € — A, + A-it +■ ( \ Z q-l -£r b ‘ tv ) 



If the delayed elastic response ( ) has associated with it a continu- 
ous spectrum of retardation times, then the factors may be represent- 
ed by a rate constant distribution function, , such that: 

b = co 

(F .4) J fl(b') ^ ~ ^ 



For this continuous distribution of retardation times the creep equation 
(F.3) becomes: 



(F.5) 



€ = A, + Alt + Ai(l- J^ftOjgT^eCb) 



Equation (F.5) now contains the Laplace transformation of the rate con- 
stant distribution function , and may be rewritten as: 



(f . 6) e - At +■ A-zt + ( 1 — 

where represents the Laplace transform of the function. 

The transform can evaluated readily in the case of 

a box distribution of rate constants. Referring to the notation used 

for the box distribution illustrated in Figure 37, 

b - 






b=o 

b=Or 



L, 

- -J 4 * 



L-jC^dLb 



b •= C, 



b*c, 



(F.7) “ 14 ) 

In equation (F.7) , constants and are the maximum and minimum rate 
constants, respectively, of the box distribution. These constants may 
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be evaluated for a given set of creep data as shown in Figure 19. 



In equations (F.l) and (F.6), the coefficient is intended to 
represent the delayed strain after elastic equilibrium is achieved, i.e., 
after infinite time. The delayed strain at time zero is zero. Therefore, 
the expression for given by equation (F.7) must be normalized 

so that: 



(F.8) XbW} = O , t - O 

< f - 9 > tWM = 1 , t « 00 



This may be accomplished by expanding the exponential terms of equation 
(F.7) in a MacLaurin series: 

tm = t [' " c ' t +<s rf - ( ^’- • - ( I -Cit + • • • )] 



(F.10) 




Ci-C, 



Equation (F.7) may now be normalized to get: 

(F.ll) = (Ci-C.Vt ^ ] 

By substitution of the time limits "t — O and *t- — 00 into equation 
(F.ll) it is apparent that the conditions of equations (F.8) and (F.9) 
are satisfied. Substitution of equation (F.ll) into equation (F.6) gives 
the analytic expression for creep strain when the retardation spectrum is 
a box distribution of rate constants. 



(F.12) 



6 = A, + A*t 



r 



+ A 




£ L|t — £ Cl t ] 
(Cx-c.vt J 
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Similar derivations may be performed whenever the retardation spectrum 
can be represented by a sum of functions whose transforms can be found. 
Unfortunately, the resulting expression is quite complex for all but the 
simplest distribution functions. For example, the wedge distribution func- 
tion, Figure 38, similar to the theoretical model of Alfrey discussed in 
Section 6, leads to the expression: 



(F- 13) 



€ = A, + A z t + A 3 



£ - k '*(k ; t— k,t-i) u ->lt ] 



In the equations for creep strain developed in this Appendix, the 
delayed elastic response, 6ol , represents a portion of the total strain. 
In equation (F.l), for example, 



(F. 14) 



= A,(\-.£' bt ) 



The value of the strain €.<£ varies from zero at "t — O to a value 
at "t- 00 . It is often convenient to normalize the delayed strain. The 
normalized form, called the delayed strain fraction, is defined by: 

(F.15) £ =• /A -j 



The delayed strain fraction varies from £ = 0 at "L~0 to £ — I at 
t — 00 , and its value at any time t > O represents the progress 

toward elastic equilibrium. For equation (F.l), the delayed strain frac- 
tion is: 



(F.16) <L - \ ~ SL 

A better insight into the nature of the delayed elastic response is 
achieved by plotting delayed strain fraction versus the logarithm of time. 
Equation (F.16) is plotted in this manner in Figure 39. The inflection 
point of this sigmoidal shaped curve is reached at a time "t — ]/u , the 
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characteristic retardation time. The delayed strain fraction at t - 
is £ — \ ^/jL 35 0,632. The slope of this curve changes very gradually 
in the range £ - 0.3 to £ = 0.8. In this region, the shape of the curve 
can be approximated by a straight line drawn with a slope corresponding to 
dt/d (X*. t) evaluated at £ = 0.5. This is shown by line segment A-A 
in Figure 39. 

For a box distribution of rate constants, the delayed strain fraction 
versus log time curve (Figure 40) is also sigmoidal in shape, with even 
less drastic slope changes in the region around £ =0.5. In order to 
accurately depict the delayed elastic response, the sigmoidal curve must 
fit the data closely at a time tj corresponding to £ =0.5. Otherwise, 
large errors in predicted strain occur over most of the time range. Fig- 
ure 41(a) illustrates a sigmoidal curve exhibiting a slope mismatch with a 
data plot. In Figure 41(b), the slopes of the data plot and the sigmoidal 
function are the same, but the responses do not occur at the same time. 

It is apparent from these figures that the error in predicted delayed 
strain fraction will be large if there is a slope and/or time mismatch at 
approximately £ =0.5. 

Plots of £ versus were made for the creep data obtained in 

this investigation. For different viscoelastic materials tested, these 
plots exhibited remarkable similarities. The delayed strain fraction 
reached a value £ = 0.5 at a time ^ 1 hour, and the slope of the 
curves at £ = 0.5 was d£/dL(U.) ^0.15. The same values were observed 
for PP, PC, NY, and PE and were independent of stress. Therefore, the 
same function may be used to characterize the delayed elastic response 
of each of these materials. A question arises as to whether or not a func- 
tion such as the box distribution leads to a suitable representation. To 
text the applicability of the box distribution, it is necessary to investi- 
gate the characteristics of the resulting sigmoidal curve for the limiting 




Delayed Strain Fraction vs. Ln Time for a Single Retardation Time 




Figure 40 

Sigmoidal Curve for a Box Distribution 





Figure 41 

Unsuitable Representations of Delayed Elastic Response 
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cases of the box distribution shown in Figure 42* This analysis is given 
below: 

If the limits of the box distribution are extended to cover many de- 
cades of time, as shown schematically in Figure 42(a), 

c z » c, 

1 ' for t >■ O 

The expression for the delayed strain fraction (from equation (F.12)) is 
given by: 



(F. 17) ^ l 

4 - 

The terms JZJ x u and 

(F . 18) 6 . - | - 



\ xr c,t - z 0 * ] 

L (c t-c.Vt J 



may be neglected so that, 




The constant is very small so that, 




The delayed strain fraction for the extended box distribution may then be 
approximated by: 

(F. 19) £ = I ~ cTt 1 Ci»C, 

This curve is illustrated in Figure 43(a) , The approximation is obvious- 
ly in error for times "t ^ •/et , but the constant C ^ 1 s large and the 
equation (F.19) is valid in the time region of interest. 

For the case shown in Figure 42(b) where the limits converge and the 
box becomes narrow, the equation for delayed strain may be approximated 
as follows: Denoting the bracketed term of equation (F*17) by 2 9 and 

expanding the exponentials, 



2 = 



(&tV (a)\ (c.fc? 

‘ U Z\ 

(C t - c,). t 



' \\ '2.'. 
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(F .20) 2 - -t -t fei-l-f 0 • V ... 

2iC^.- c 0 3!CCi' c i) 

'Ai, 

Th? torn of tills expansion is represented by: 

C/" — Cb^ , 

The value of this term for the case of converging limits C C ^ 

may be determined by considering C 9 and C n to differ by a small quantity £ 






r ~ r 
^2 1 



then 



C, 



'Vv 



(C 1 * ; 



/TL 



Using a binomial expansion and truncating after the second tern, 

i 



M 



1 



' v«V/ ^ M * 

= (c.4 e) - c* + ^et| j 4 



then. 



cS-c* 



mCC ( 



+* ~ f 



^ r Rtf, 



r M- l 

_ C, 



/V\ 



.-•‘Uf. *c-0 ^ ! £ (/n-ut 



therefore, 



(F.21) 



r* r A-v 

* •' h 



l v ) 

c* c > a\{ (Y* ~ c 



c 



"V-/ 



(aa~!) 



Substituting equation (F.21) for the m^tem of equation (F.20) aixd 
evaluating the terms for /XL = 1,2,3 the expansion becomes: 



(F.22) 7 = 



C.t 

I - “',7 



* < 7 Al — J 

r,‘t __ C, £*•’• 

Tl 



„-c,t 

Eouation (F.22) represents the expansion of J2. in a MacLaurin 

series, so the desired approximation has bean found: 

-< ,t 



(F . 23) 



f ~ l - Ji 



or 



<f * I “ * 



- c.. 



for C, 
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Equation (F.23) would result from a n single spike" distribution, corres- 
ponding to a single characteristic retardation time. The analysis veri- 
fies the expected result that equation (F.17) should reduce to equation 
(F.23) as the width of the box distribution becomes vanishingly small. 

This approximation is sketched in Figure 43(b). 

In order to be used for an accurate representation of the data in this 
investigation, the slope of the sigmoidal plot of the box distribution must 
take on the value: 



(F.24) 



dz 



£=0.5 
t - I 



o. \s 



By the approximation for the limiting case of C^, given by equa- 

tion (F.19), 

t = \ ~ 

for ^ O.S } "t — l 

C x = Z 

sLE _ -L . - dL.t 

dJ. _ _L 



(F. 25) 



ole 

d(Uit) 



£ = 0*5 
\ 



0.5 



For the case of narrowing box width, (C^ — *■ C^), given by equation 

(F.23) : 

t- I - JL~ 




d. t 




for £ = 0.5 , t - I 

c^^o.7 
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c x »> C, 



9ft>) 




Cl 



b — 






Figure 42 

Limiting Cases of the Box Distribution 




Figure 43 

Approximations of the Response Function for the 
Limiting Cases of the Box Distribution 
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(F. 26) 





It may be demonstrated by numerical methods that equation (F.26) represents 



(F.24) can not be met, since the observed slope of 0.15 does not fall in 
the range of 0.35 to 0.5. Therefore the box distribution is an unsatis- 
factory choice for the retardation spectrum of the materials tested in 
this experiment. 

A similar analysis of the wedge distribution (Figure 38) shows a 



indicating that this rate constant distribution is also unsatisfactory. 

Since the limiting case of the box distribution given by equation (F.23) 
corresponds to the response of a single Voigt Element (Figure 5), the four 
parameter creep model (Equation (F.l)) gives a poor representation of de- 
layed elastic strain. 

As discussed in Section 6, good correlations can be obtained using 
a model with several discreet retardation times. This 'Wilti-spike distri- 
bution" produces a plot of delayed strain fraction versus log time with 
inflections at each retardation time. These retardation times may be chosen 
so that the response function approximates the creep data quite closely. 
Figure 44 illustrates the response function resulting from a three-spike 
distribution superimposed on a typical plot of creep data. The response 
function from a box distribution is also shown. The lower limit, C^, of 
the box is equal to the slow rate constant b>j of the three-spike distri- 
bution, and the upper limit equals the fast rate constant . 



the minimum value of the slope for any choice of the 

constants and C^ 9 provided The requirements of equation 





0.4 for all choices of and K^, 



115 



These values for the box distribution function were chosen to illustrate 



the dominant role of the fast rate constants in governing the delayed 
elastic response. The constants and C ^ could be changed to shift the 
sigmoidal curve toward longer times, but the slope cannot be adjusted to 
produce a good fit with the data. 

Monomodal distribution functions characterized by a single maximum 
can always be approximated by a sum of box- type distributions as shown in 
Figure 45. The limits of each succeeding box fall within the range of the 
first box. Since the slope of the delayed strain fraction versus log 
time plot resulting from such a distribution is restricted in accordance 
with equations (F.25) and (F.26), it is apparent that the retardation 
spectrum of the plastics tested is not monomodal. 

Bimodal (or polymodal) spectra produce a series of inflections in the 
plot of £_ vs. , which results in a lower overall slope, i.e., the 

delayed elastic response occurs over many decades of time. The three- 
spike distribution of Figure 44 is an example of a polymodal retardation 
spectrum. 

The delayed elastic response of the materials tested covered about 
five decades of time. This response can be characterized by a polymodal 
retardation spectrum such as illustrated in Figure 46, where the peaks 
imply the presence of several different molecular mechanisms whose time 
constants differ by orders of magnitude. This model is in accord with 
Alfrey f s concept (See Section 6) of local bond straightening and long- 
range uncoiling of the polymer chains. 



116 



00 




* 

u 

3 

o 

X 



of 

£ 

£ 

c 



117 



Figure 44 

Delayed Elastic Response Functions Fitted to Creep Data 




Figure 45 

Monomodal Retardation Spectrum Approximated By Box Distributions 




Figure 46 

A Polymodal Distribution Function 
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APPENDIX G 



EQUIPMENT 

The creep machines used in this investigation apply constant dead- 
weight loading through calibrated lever arms. (Figure 47). Eleven identi- 
cal units are installed in one test stand. Photographs of the laboratory 
and apparatus are shown in Figures 48 through 50. Specimens were gripped 
by steel pins placed through the 5/8 inch diameter hole in each end, and 
uniaxial loading was maintained through a series of universal joints 
above and below the grips. Turnbuclcles and an adjustable counterbalance 
weight permitted final balancing of the lever arm prior to load applica- 
tion. 

Studies of viscoelastic behavior of materials require large amounts 
of data obtained over long time periods of testing. Consequently, facili- 
ties for automatic recording of creep and relaxation data are highly de- 
sirable. Design of such instrumentation was started in conjunction with 
this project and is presently being installed. Special clip-on strain 
gage fixtures were fabricated to permit measurement of the large strains 
encountered in mechanical testing of viscoelastic materials. These clip 
gauges, shown in Figures 51 and 52 consist of phosphor bronze strips with 
BLH SR-4 Type A-5-1 strain gauges mounted on each side. Tests on phosphor 
bronze indicate that this material has remarkable linearity over the en- 
tire elastic range with no observable hysteresis effects. The strips are 
connected to aluminum legs which make a knife-edge contact with the speci- 
men, This design has two useful features: 

(1) The legs are detachable so that different leg lengths may 
be used depending on the strain sensitivity desired. 

(2) The mode of attachment to the specimen is such that this 
electromechanical device may be used for specimens of various 
widths and thicknesses. 
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The four strain gauges are connected in a Wheatstone bridge arrange- 
ment shown in Figure 53. A potentiometer in the circuit permits a null 
adjustment to be made. This bridge circuitry has several advantages: [26] 

(1) Automatic moment compensation (in the event there are devia- 
tions from axial loading. 

(2) Signal strength increased by a factor of four, (over that 
of a single gauge mounted on the specimen) 

Terminal boards have been set up on the eleven-unit creep machine to 
receive the leads from the extensometer circuit for each test specimen. 
Bridge circuits, each with its own balancing potentiometer, are placed 
in parallel across a D.C. power supply, with leads from each bridge enter- 
ing master control panels. From these panels, the output signal from 
each extensometer may be fed into a multi-channel recorder, permitting 
testing and graphical data output for up to eleven specimens simultaneous- 
ly. A stress analysis of the clip gauges, using linear elastic theory, 
is given in Figures 54 and 55. 

It is desirable to have the output of electromechanical extensometers 
converted directly into digitalized form because of the usefulness of the 
digital computer in data analysis. Preliminary design of such instru- 
mentation has begun at NPGS, Monterey. 
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counter balanca. 




Figure 47 

Lever Arm Loading Assembly 




Figure 48 

dpecimen Grips and Clamps 
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Figure 49 

Eleven f nit reep Testing Apparatus 




Figure 50 



Tripoli tis Extensometer Attached to Specimen 
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Clip Gauge Extensometer 
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Clip Gauge Extensometer Mounted on Specimen 
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Figure 53 

Bridge Circuitry for Extensometer 
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Figure 54 

Dimensioned Sketch for Clip Gauge Stress Analysis 
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CLIP GAUGE STRESS ANALYSIS 



Assume no bending in aluminum legs 

To find E in gauges for a given specimen extension, S 

e = 



pe-T 



from the flexure formula, £ — 'i/e 

P“ , <f=44 



^ (lh) 

for small angles 0 , XitC* 

£ ■= ^VzLh 

for d = .Olo", b«2", L-l" 



s 

2b 



6 = .00125 s 



for a specimen extension of $ — | M 

6 - l250/<.»ncb«/incb 

or about one -ha If the maximum strain 
capability of the BLH type a- 5-1 gages 



maximum stress in the bronze strip: 



CT- Ee , 



€~.OOtT.S 

E- * loops’* 



<y - 2.0 poo p*\ 

or only a small fraction of the yield 
strength of phosphor bronze ( Oy = 112,000 psi) 

at an extension of S“ I " , error in the strain 
due to the small angle assumption is: 

%h= Al «-tS , W =.2474 

• Error- .. 



Figure 55 

Theoretical Stress Analysis of Clip Gauges 
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